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FOREWORD 


This  report  deals  with  the  prediction  of  explosion  injury  to  fish  with 
swim  bladders  and  is  part  of  a  continuing  study  of  the  effects  of 
underwater  explosions  on  marine  life.   Swim  bladder  fish  are  parti- 
cularly vulnerable  to  explosions,  and  this  group  includes  the 
majority  of  fish  with  sports  and  commercial  value.   This  study  will 
result  in  an  improved  capability  to  predict  such  effects,  and  will  be 
useful  in  connection  with  the  testing  of  new  explosives  and  of  warheads 
at  sea. 

This  study  is  part  of  the  ordnance  pollution  abatement  program  of  the 
Naval  Sea  Systems  Command  and  was  supported  by  SEA  TASK  SSL  55001/ 
19373. 

The  author  is  indebted  to  Ermine  A.  Christian  and  George  A.  Young  for 
many  valuable  suggestions  during  the  course  of  this  work. 

JULIUS  W.  ENIG 
By  direction 
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DTNEG 


Duration  of  the  negative  pressure  phase  (Figure  3.1.1) 
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NSWC/WOL  TR  76-155 

PMAX     Shockwave  peak  overpressure  referenced  to  initial 
ambient  pressure,  p.,  but  generally  referred  to  as 
"Shockwave  peak  pressure" 

0         Shockwave  decay  constant  (Figure  3.1.1) 

PCORR     Coefficient  used  to  adjust  Shockwave  peak  overpressure 
for  batch-to-batch  variations  in  explosive  charges 
(Equation  3.1.1) 
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p„   (t)   Shockwave  overpressure  calculated  by  superposition 
of  direct  and  surface-reflected  waves, 

PSUM(t)  =  PD(t)  +  PR(t) 
PATM     Atmospheric  pressure 

y        Depth  of  fish 

y.        Depth  of  top  of  cavitation  (on  surface-reflected  ray) 

y        Depth  of  cavitation  closure  point  (on  surface-reflected 
ray) 


R         Slant  range  from  charge  to  cavitation  closure  point 

PMAX      Shockwave  peak  overpressure  at  cavitation  closure  point 

6        Shockwave  decay  constant  at  cavitation  closure  point 

R„        Slant  range  from  charge  to  point  of  reflection  at 
water  surface 


DOB 


Depth  of  burst 
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DYNAMICAL  MODEL  FOR  EXPLOSION  INJURY  TO  FISH 

1    INTRODUCTION  AND  SUMMARY 

This  report  describes  a  method  for  calculating  the  kill  probability 
for  bladder  fish  subjected  to  an  underwater  explosion.   The  method 
consists  of  an  approximate  calculation  for  the  extreme  values  of  com- 
pression and  extension  of  the  fishes'  gas-filled  swim  bladder  in 
response  to  the  explosion  pressure  wave.   The  calculations  are  made 
for  the  damped  radial  oscillations  of  a  spherical  air  bubble  in  water. 
The  kill  probability  is  then  calculated  as  an  experimentally  determined 
function  of  the  ratio  of  maximum  to  minimum  radius  during  the  oscilla- 
tory response.   For  each  species  of  fish  the  effective  swim  bladder 
radius  needed  for  the  calculations  is  determined  by  trial  and  error 
thru  correlating  the  observed  injuries  with  the  calculated  ratio  of 
maximum  to  minimum  radius. 

The  method  is  used  to  describe  dissection  results  from  1500  Spot 
and  White  Perch  caged  at  depths  from  5  to  100  feet  and  subjected  to 
pentolite  explosions  of  from  1  to  70  lbs  submerged  at  depths  ranging 
from  from  5  to  70  feet.1'2 

These  results  are  then  used  to  calculate  several  examples  of 
kill-probability  and  fish-kill  distributions  for  a  particular 
explosion  geometry  and  several  species  and  sizes  of  fish.   Details 
of  the  computations  are  given  in  the  appendices.   Appendix  C  gives 
an  approximate  method  for  calculating  the  pressure  signature  from  an 
underwater  explosion  subject  to  surface  effects  including  cavitation. 


1.  Gaspin,  Joel  B.,  "Experimental  Investigations  of  the  Effects  of  Underwater 
Explosions  on  Swim  Bladder  Fish,  I:   1973  Chesapeake  Bay  Tests",  NSWC/WOL  TR  75-58, 
1975. 

2.  Gaspin,  J.  B.,  Wiley,  M.  L.,  and  Peters,  G.  B.,  "Experimental  Investigations  of 
the  Effects  of  Underwater  Explosions  on  Swim  Bladder  Fish,  II:  1975  Chesapeake  Bay 
Tests",  NSWC/WOL  TR  76-61,  19  76. 

13 


NSWC/WOL  TR  76-155 


DYNAMICAL  MODEL 


2.1  THE  SWIM  BLADDER.   Figure  2.1.1  shows  a  life-size  tracing  of 
a  Spot  (Leiostonms  xanthurus  Lacepede)  which  depicts  the  swim  bladder 
and  the  surrounding  internal  organs.   Figure  2.1.2  shows  a  swim  bladder 
removed  from  a  Spot.   In  our  work  we  depicted  the  motion  of  this  organ 
by  that  calculated  for  a  spherical  bubble  of  air.   The  size  of  this 
"equivalent  spherical  bubble"  was  first  estimated  by  equating  its 
volume  to  that  of  the  fishes'  swim  bladder  —  then  later  we  adjusted 
the  size  of  the  equivalent  bubble  to  optimize  the  correlation  between 
the  calculated  motion  and  the  observed  injuries. 

Figure  2.1.3  shows  a  life-size  tracing  of  a  White  Perch 
(Mo rone  americana  Gmelin) ,  the  other  species  of  fish  used  in  the  work 
reported  here.   The  inset  shows  the  shape  of  the  forward  portion  of 
the  swim  bladder  in-place  in  the  fish. 

2.2  PRELIMINARY  CALCULATIONS  OF  BLADDER  OSCILLATION.   We  began 
this  study  by  making  order  of  magnitude  calculations  which  approxi- 
mated the  explosion  pressure  field  with  a  step  increase  followed  by  a 
step  decrease  in  pressure.   The  fishes'  swim  bladder  was  approximated 
by  a  spherical  bubble  of  air  in  an  infinite  body  of  water  —  and,  we 
calculated  its  undamped  oscillatory  response.   These  order  of  magnitude 
calculations  showed  that  the  oscillatory  response  of  the  swim  bladder 
was  a  likely  source  of  the  fishes'  injuries.   They  also  pointed  out  a 
strong  resonance  which  occurs  when  surface  cut-off  (the  arrival  of  the 
rarefaction  wave  reflected  from  the  water  surface)  happens  at  the 
instant  of  maximum  bladder  compression.   This  work  was  described  in  an 
internal  Center  report  which  is  reproduced  as  Appendix  A  of  this  report. 

2.3  CURRENT  METHOD  FOR  CALCULATING  BLADDER  OSCILLATION.   The  method 
developed  for  calculation  of  the  oscillatory  response  to  the  changing 
pressure  field  generated  by  an  underwater  explosion  is  described  in 
detail  in  Appendix  B.   Here,  we  attempt  to  give  the  reader  a  feel  for 
this  computation  and  its  approximate  nature.   Figure  2.3.1  shows  the 
idealized  pressure  signature  measured  on  an  underwater  explosion  test 
together  with  the  calculated  pressure  inside  a  fish's  swim  bladder 

14 
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FIG.   2.1.2     SWIM   BLADDER    REMOVED    FROM  SPOT 


16 


NSWC/WOL  TR  76-155 


I 
o 


to 

O 


w 

UJ 

u. 


I 
o 

UJ 

a. 

UJ 

X 

u. 
O 

O 

2 

o 
< 

a. 

> 
a 
o 
a 

CO 

ei 


UJ 


17 


NSWC/WOL  TR  76-155 


Nt 

LU 

> 
< 
§ 
LU 
CC 
D 

LU 

Ol 

CC 
a. 

z                            ^S 
o                          ^^ 

CO                                 f 

o                   V 

o 

LU 

s 

H 

LL 

o 

cc 

Q- 

X                            r 

LU  -^....^^^ 

o 

z 
o 

o 

z 

D 

LL 

LU 
Q 
Q 

< 

_i 
m 

z 

o 

< 
CO 
< 
LU 

N 

CO 

LU 
DC 

D 
CO 
CO 
LU 
CC 
0. 

N                             LU 
/I                           CO 

o 

TIME 
5SURE  an: 

\_ 

/     J                             ^~ 

y           lu 

IS                            Q 

^--^/                          < 

— — "              /                          _i 

/                           m 

o 

LU 

CC 

OL 

DC 
LU 

Q 

o 

xo 

_l 

CO 
CO 

0    j 

O 

u. 

3anss3dd 


18 


NSWC/WOL  TR  76-155 

as  it  responds  to  this  applied  pressure.   Below  the  pressure  traces  are 
corresponding  sizes  calculated  for  the  spherical  air  bubble  used  to 
represent  the  fish's  swim  bladder. 

This  fluid  motion  can  be  likened  to  a  mass  attached  between 
two  springs,  the  end  of  one  being  fixed,  the  end  of  the  other  movable. 


'ff^MSu 


'^-vJlQJL/-^^ 


The  first  spring  is  analogous  to  the  internal  gas  pressure;  the 
second  to  the  changing  external  water  pressure.   Let  us  imagine  the 
mass  resting  on  a  frictionless  block  of  ice.   It  is  at  rest  in  a 
position  of  equilibrium.   If  at  t=0  we  quickly  displace  the  movable 
spring  to  a  new  position,  the  mass  will  oscillate  about  a  new  position 
of  equilibrium  determined  by  the  displacement  of  the  end  of  the  spring, 
The  mass  will  oscillate  about  this  new  equilibrium  position  with  an 
amplitude  equal  to  the  displacement  of  this  equilibrium  position  from 
the  initial  at-rest  position  of  the  mass. 

This  initial  at-rest  position  now  becomes  an  extremum  of  the 
oscillatory  motion.   Since,  at  each  extremum  of  the  motion  the  mass 
is  instantaneously  at  rest,  a  subsequent  jump  displacement  of  the  end 
of  the  spring  occurring  at  an  extremum  simply  shifts  the  motion  to  an 
oscillation  about  a  new  equilibrium  position  -  and  the  new  amplitude 
is  just  the  distance  of  the  mass  at  the  time  of  the  jump  from  this 
new  equilibrium  position. 

Likewise,  for  the  bubble,  step  changes  in  the  outside  water 
pressure  occurring  at  half-period  intervals  simply  change  the  equili- 
brium pressure  (equal  to  the  outside  pressure)  of  the  oscillating 
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bubble  flow  -  and  the  new  amplitude  is  just  the  instantaneous 
excursion  of  the  bubble  pressure  from  the  changed  outside  pressure. 
We  made  use  of  this  property  of  the  motion  to  calculate  the  motion  as 
the  response  to  a  series  of  pressure  steps  as  shown  in  Figure  2.3.2. 
(Appendix  B  describes  the  way  we  estimated  successive  half-periods 
in  order  to  construct  the  approximating  step  wave.) 

The  transition  from  the  positive  to  the  negative  pressure 
phase,  i.e.,  surface  cut-off,  must  be  treated  as  a  special  case,  since 
in  general,  this  will  not  occur  at  an  extremum  of  the  motion.   If 
cut-off  occurs  at  a  compression,  this  compression  is  then  the  compres- 
sion of  the  next  half  cycle  also  -  and  the  subsequent  expansion  is 
the  greatest  possible.   If  cut-off  occurs  at  an  expansion  and  the  new 
outside  pressure  is  lower  than  that  inside,  this  expansion  becomes 
the  compression  of  the  next  half  cycle,  there  occurs  a  180  degree 
phase  shift,  and  the  subsequent  expansion  is  the  smallest  possible. 
In  general,  cut-off  occurs  somewhere  in  between  and  there  is  an 
intermediate  phase  shift  and  subsequent  expansion  or  compression. 
In  such  cases  we  calculated  this  next  half  cycle  at  lowered  pressure 
as  starting  at  the  time  of  the  previous  extremum  but  at  a  ficticious 
amplitude  interpolated  as  a  function  of  the  time  at  which  cut-off 
occurred.   (The  details  of  this  procedure  are  described  in  Appendix  B.) 

Finally,  at  the  end  of  the  negative  phase,  the  return  to 
ambient  pressure  can  cut  short  or  even  prevent  the  final  expansion. 
This  was  handled  in  an  analogous  manner  to  surface  cut-off  and  is 
also  described  in  Appendix  B. 

To  approximate  the  dissipation  of  energy  which  in  nature  must 
occur  during  the  oscillation,  we  arbitrarily  extracted  a  fraction  — 
30%  —  of  the  remaining  oscillation  energy  at  each  half  cycle  of  the 
motion. 
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3.   CORRELATIONS  WITH  EXPERIMENTAL  DATA 

3.1  SUMMARY  OF  TEST  DATA  AND  COMPUTATION  RESULTS.   Tables  3.1.1 
thru  3.1.4  summarize  the  input  and  output  parameters  of  the  computa- 
tions of  this  report.  These  tables  include  all  the  test  data  for  Spot 
and  White  Perch  obtained  on  the  1973  and  1975  test  series  in  the 
Chesapeake  Bay.   For  the  1973  test  data  (the  "500"  Shot  Numbers)  the 
fishes'  horizontal  ranges  and  depth  which  are  listed  in  the  tables 
are  the  preset  or  intended  values.1   On  the  1975  tests,  shots  783 
thru  786,  water  currents  caused  significant  deviations  in  horizontal 
range  and  depth.   For  these  shots  the  preset  values2  have  been 
corrected  as  described  in  Appendix  D.   (Tables  D-2  thru  D-7  list 
corrected  values  for  charge  depth  and  fish  and  pressure  gage  locations 
for  all  shots  of  the  1975  test  series.) 

The  ambient  pressure  (Tables  3.1.1  and  3.1.2)  is  the  total 
hydrostatic  head  at  fish  depth  calculated  for  fresh  water,  density  = 
1000  kg/m3,  and  nominal  atmospheric  pressure,  1.013  x  10 5  pascals. 
The  vapor  pressure  was  calculated  from  the  measured  water  temperature 
(1975  tests)  or  estimated  water  temperature  (1973  tests) . 

The  peak  pressure  PMAX  in  pascals  and  time  constant  0  in 
milliseconds  are  calculated  from  the  similitude  relations  for 
pentolite 


PMAX  =  PCORR  x  5.  62 


/T7l/3\  1.14 
x  107  ^L_ J  (3.1.1) 


'1/3  fcfr) 


0  =  0.084  Wx/J  l-^\    °-23  (3.1.2) 


where  W  is  the  explosive  mass  in  kilograms,  R  is  the  slant  range  from 
the  charge  in  meters,  and  PCORR  is  an  empirical  constant  that  was  used 
to  adjust  Equation  3.1.1  for  batch-to-batch  variations  in  the  explo- 
sive charges.   For  the  1973  test  series  PCORR  =  1.10  —  for  the 
1975  tests  PCORR  =  0.95. 

(Text  continued  on  naae  32) 
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The  positive  duration  TPOS  (see  sketch  Figure  3.1.1)  was 
calculated  as  described  in  Appendix  C  for  the  1973  tests;  for  the 
1975  tests  measured  values  were  used.   The  negative  pressure  PNEG 
and  its  duration  DTNEG  were  measured  from  the  pressure  records. 

The  last  five  columns,  Tables  3.1.1  and  3.1.2,  list  calculated 
parameters  for  the  oscillatory  response  of  the  swim  bladder.   The 
number  of  positive  cycles  is  the  positive  duration  measured  in 
complete  cycles  of  oscillation  starting  from  the  initial  "at  rest" 
expansion.   The  number  of  negative  cycles  locates  the  final  return  to 
ambient  pressure  (end  of  negative  phase  —  see  Figure  3.1.1)  in  terms 
of  the  bladder  oscillation.   "Zero"  is  taken  as  the  last  expansion 
occurring  before  surface  cut-off.   See,  e.g.,  Figure  2.3.1.   Thus, 
unless  the  number  of  negative  cycles  is  greater  than  0.5,  no  final 
expansion  takes  place  during  the  negative  phase.   The  number  of 
negative  cycles,  (At/T)   ,  is  a  measure  of  the  contribution  of  the 
negative  phase  to  the  final  bladder  expansion.   The  relationship  is 
roughly  as  follows: 


(At/T)    <  0.5  No  contribution 


0.5  <  (At/T)„_,  <  1.0  Partial  contribution 

NP 


1.0  £  (At/T)  Full  contribution 


The  reader  is  referred  to  Appendix  B  for  a  more  complete  description, 

The  bladder  oscillation  parameters  X,  Y,  Z  are  defined  as 
follows : 


X  =  -100  In  AMIN/A.  (3.1.3) 


Y  =  10  0  In  AMAX/A.  (3.1.4) 

Z  =  X  +  Y  =  100  In  AMAX/AMIN  (3.1.5) 
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where  In  is  the  natural  logarithm,  Aj^  is  the  initial  at  rest  bladder 
radius,  and  AMIN  and  AMAX  are  the  smallest  and  largest  radii  during 
the  oscillatory  response.   In  the  cases  studied  so  far,  the  first 
compression  has  always  been  AMIN,  while  AMAX  generally  occurs  during 
or  following  the  negative  phase.   We  will  use  the  bladder  oscillation 
parameters  X,  Y,  and  Z  in  Sections  3.2  and  3.3  to  discuss  and  correlate 
the  observed  injuries  to  fish  specimens  on  our  tests. 

The  bladder  radius,  A.  —  Tables  3.1.3  and  3.1.4,  column  7  — 
is  the  radius  for  a  sphere  of  the  same  volume  as  the  fishes'  bladder 
at  fish  depth  for  a  fish  which  has  not  acclimated,  i.e.,  for  a  fish 
at  equilibrium  with  one  atmosphere  pressure  which  is  isothermally 
compressed  to  hydrostatic  pressure  at  fish  depth.*  Bladder  radius  A. 
for  non-acclimated  fish  is  given  by 

A.  =  (A.)    A^Oy_x_loA  V3  (3.1.6)** 

1  °   V    pi     / 


where  p.  is  the  pressure  at  fish  depth  in  pascals  and  (A.)  ,  the 
radius  at  one  atmosphere  pressure,  is  given  by 

(A±)  /L  -  0.033  (Spot)  (3.1.7) 


(A.)    /L   =    0.055  (White    Perch)  (3.1.8) 


where   L    is    the    fork    length   of   the    fish.      The   constants    0.0  33   and 
0.055   were    selected    for   each    species   by   trial   and   error   so   as    to 
optimize   the   correlation   between   experimentally   observed   injuries    to 
the    test    specimens    and    the    calculated   bladder   oscillation    parameter 
Z     (e.g.,    Figure    3.2.1). 


*This    is   how  it  was   done   on   the   explosion   tests.      It  would  have  been  prohibitively 
expensive    to   have   let    the    fish   acclimate    to    testing  depth,    because  of   the   time 
required.      This   problem  is   discussed  by  Gaspin,   Wiley,    Peters,    1976. 

**  For   the  calculations   of   this   report    the   exponent,    l/3y,    where  Y  -   1.4   —   for 
adiabatic  compression  —  was    inadvertently   used.      Since   the  computations   reported 
here  should  be  quite   insensitive   to    this   error,    they  were  not   redone. 
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The  number  of  fish  of  a  given  average  size  at  each  test  specimen 
location  and  summaries  of  the  injuries  evaluated  upon  post-shot 
disection  are  listed  in  the  last  5  columns  of  Tables  3.1.3  and  3.1.4. 

Thus,  e.g.,  on  Shot  531  (Table  3-1.3)  at  the  76.2  meter  range 
position  there  were  10  Spot,  all  10  of  which  received  injuries  of 
level  1  or  greater,  7  of  which  received  injuries  of  level  2  or  greater, 
5  of  which  received  injuries  of  level  3,  and  none  of  which  received 
level  4  injuries.   (There  were  no  injuries  of  level  5  observed  on 
either  the  1973  or  1975  tests.) 

Table  3.1.5  is  a  condensed  version  of  the  code  used  to  evaluate 
the  explosion  damage  to  the  fish  upon  post-shot  dissection.3 
Successive  injury  levels  are  of  increasing  severity  and  each  fish  is 
classified  to  fall  into  one  and  only  one  level  of  injury.   Figure 
3.1.2  depicts  this  fundamental  design  of  the  injury  code.   Thus,  the 
fundamental  physical  problem  of  this  report  is  to  estimate  test  condi- 
tions for  the  transitions  from  one  injury  level  to  the  next.   We  will 
do  this  by  estimating  the  probability  of  events  such  as  "occurrence  of 
injury  level  3  or  greater"  as  a  function  of  the  calculated  damage 
parameter  Z  (Equation  3.1.5).   It  is  convenient  —  but  not  fundamental 
to  the  analysis  --  that  such  "cumulative"  injuries  are  also  what  is 
of  practical  interest  in  explosion  testing.   For  example,  the  event 
"injury  of  level  3  or  greater"  has  been  found  to  correlate  with  the 
observed  fish-kill*  on  underwater  explosion  tests.  And,  the  event 
"injury  of  level  2  or  greater"  has  been  considered  to  give  the  fish 
"little  chance  to  survive  predation  .  " 


3.   The  injury  evaluation  code  is  due  to:   Hubbs ,  C.  L. ,  Schultz,  E.  P.,  and 
Wisner,  R. ,  Unpublished  preliminary  report  on  "Investigation  of  Effects  on  Caged 
Fishes  from  Underwater  Nitro-Carbo-Nitrate  Explosions,"  U.  of  California,  Scripps 
Institute  of  Oceanography,  1960.   The  complete  code  is  also  listed  in  reference  1. 

*  The  term  "observed  fish-kill"  refers   to  the  dead  and  dying  fish  found  on  the 
surface  and  bottom  following  an  underwater  explosion. 
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Injury 
Level 


TABLE  3.1.5 


CODE  FOR  EVALUATING  VISIBLE  FISH  DAMAGE 


No  damage 

Light  hemorrhaging  in  tissues  covering  kidney 

Light  hemorrhaging  throughout  body  cavity,  some  kidney 
damage,  but  gas  bladder  intact 

Severe  hemorrhaging  throughout  body  cavity,  gross  kidney 
damage,  and  gas  bladder  burst 

Partial  break-thru  of  body  wall,  bleeding  about  anus 

Ruptured  body  cavity,  internal  organs  scrambled  or  lost 
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However,  should  the  reader  require  for  example  the  probability, 
P (level  2  injury),  it  is  just  the  difference,   P (level  2  or  greater) - 
P (level  3  or  greater). 

3.2   OBSERVED  INJURIES  TO  SPOT  VS  CALCULATED  OSCILLATION  PARAMETER 
Figure  3.2.1  shows  the  experimentally  observed  injuries  (as 
shown  in  each  row  of  Table  3.1.3)  to  Spot   on  the  1973  and  1975  tests 
plotted  as  a  function  of  Z  =  100  In  AMAX/AMIN.   There  are  three  plots, 
one  for  each  of  the  injury  levels  1,  2,  and  3.   The  plotted  points 
represent  the  percent  of  fish  of  a  given  size  at  a  given  specimen 
location  receiving  injuries  of  the  indicated  level  or  greater. 

To  investigate  the  functional  dependence  of  the  observed 
injury  on  the  calculated  oscillation  parameter  Z  we  constructed 
Table  3.2.1  and  Figure  3.2.2.   Table  3.2.1  was  constructed  from 
Table  3.1.3  by  reordering  (sorting)  the  entries  in  order  of  increas- 
ing value  of  Z  and  then  separating  the  new  table  into  groups  repre- 
senting approximately  100  fish.   The  total  observed  injuries  for  each 
group  was  then  summarized  by  the  successive  row  entries  in  Table  3.2.1. 
Thus,  for  the  first  group  with  a  mid-range  value  for  Z=54;  41%  of  the 
fish  received  injuries  of  level  1  or  greater;  11%,  level  2  or  greater; 
and  5%,  level  3  or  greater. 

Figure  3.2.2  shows  plots  of  the  averaged  injury  data  tabulated 
in  Table  3.2.1.   The  solid  curves  are  drawn  by  eye  through  the  data 
points.   The  experimental  data  for  level  3  injuries  to  Spot  exhibit 
less  scatter  than  that  for  levels  1  or  2 ;  and  also  less  than  for  all 
levels  of  injury  to  White  Perch  (Section  3.3),   Whether  or  not  this 
has  any  physical  significance  is  not  known. 

Also  shown,  as  a  dashed  line  in  Figure  3.2.2,  is  the  equation 

(3.2.1) 


1  +  EXP [-0.055(Z-125) ] 


where  p  is  the  probability  of  observing  injury  of  level  3  or  greater, 
Equation  3.2.1  is  the  Cumulative  Logistic  Probability  Function  with 


38 


MSWC/WOL  TR  76-155 


100% 


50    - 


0 
100% 


o 

-]0— o— 

a 
a 
a               a         o 

o                       a 

— po-oo — o ooooo o □ 

- 

o 

□                     DO 

- 

D 

o 

o 

o° 

a 

INJURY  LEVEL  [T) 

o-o 

1 

— o o — bo— Q 

i                        i 

a 
in 

3 


I 

LL 


C3 
< 


O 

cc 


50 


-ooooo o o 


□  CD      O 


INJURY  LEVEL  [2] 

J 1 


100% 


50 


□  00         □ 


INJURY  LEVEL  [3] 


50  100 


150 


200  250 


FIG.  3.2.1    OBSERVED  INJURIES  TO  SPOT  AS  A  FUNCTION  OF  Z 
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FIG.  3.2.2  OBSERVED  CUMULATIVE  INJURY  TO  SPOT  AS  A  FUNCTION  OF  Z 
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the  parameter  values,  0.055  and  125,  fitted  to  the  experimental  data. 
The  parameter  value  0.055  was  obtained  by  considering  all  of  the  Spot 
and  White  Perch  data,  levels  1,  2  and  3  from  the  two  test  series. 
The  parameter  value  125  --  the  value  of  Z  corresponding  to  50% 
injury  probability  --  is  the  maximum  likelihood  fit  value  to  the 
level  3  injury  data  for  Spot. 

Equation  3.2.1  is  considered  to  be  a  sufficiently  precise 
representation  of  the  data  for  most  applications,  such  as  those 
described  in  Section  4  of  this  report. 


3.3   OBSERVED  INJURIES  TO  WHITE  PERCH  VS  CALCULATED  OSCILLATION 
PARAMETER 

Figure  3.3.1  shows  the  experimentally  observed  injuries  to 
White  Perch  on  the  1973  and  1975  tests  (as  shown  in  each  row  of 
Table  3.1.4)  plotted  as  a  function  of  Z .   To  investigate  the 
functional  dependencies  we  constructed  Table  3.3.1  and  Figure  3.3.2 
from  the  data  presented  in  Table  3.1.4.   The  plot  shown  in  Figure 
3.3.2  for  level  3  injuries  exhibits  excessive  scatter.   Perhaps 
these  level  3  injuries  to  White  Perch  are  more  closely  related  to 
some  other  parameter  of  the  bladder  oscillation,  but  it  did  not 
appear  worthwhile  to  investigate  that  possibility. 

Figure  3.3.3  shows  the  averaged  injury  data  for  both  Spot  and 
White  Perch  (Tables  3.2.1  and  3.3.1)  plotted  as  a  function  of  Z.  The 
curves  were  drawn  by  eye  through  the  combined  data.   These  plots 
support  the  viewpoint  of  this  report  that  the  observed  injuries  to 
both  species  of  fish  --  Spot  and  White  Perch  --  are  described  by  the 
same  damage  functions.   (The  coefficients  —  Equations  3.1.7  and 
3.1.8  --  for  the  effective  bladder  radii  used  in  the  dynamical 
response  calculations  are  different  for  the  two  species  of  fish, 
however. ) 


42 


NSWC/WOL  TR  76-155 


100% 


V) 


% 


o 
cc 

Q. 


50 


0 
100% 


-CDtpXEDO- 


50 


0 
100% 


50 


INJURY  LEVEL  (T) 

J I 


— i *ppi«»-rrTW»fWwTTT>-«-T-w-M-rTrwTr>on o JiiiiMiitn o — O 

□  O 


□ 

a 


INJURY  LEVEL 


OO-i o 

j— oo oap — nm oo 

o  a 

a 

a 

□ 

V 

a 

D 

a 

ao                                                                   — 

a 

G 

□ 

o                0 

a 

°        8 

o                 o 

a 

INJURY  LEVEL  @) 

a 

n-TTTHTJ ' 

50 


100  150 

z 


200  250 


FIG.  3.3.1   OBSERVED  INJURIES  TO  WHITE  PERCH  AS  A  FUNCTION  OF  Z 
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3.4   COMPARISON  OF  SHOCKWAVE  IMPULSE  AND  CALCULATED  BLADDER 
OSCILLATION  AS  DAMAGE  PARAMETERS 

Yelverton  et  al4  at  the  Lovelace  Foundation  working  in  a  test 
pond  with  eight  species  of  bladder  fish  at  depths  down  to  10  feet 
obtained  fish-kill  results  describable  by  an  Impulse  Damage  Parameter, 

Z   =  80  In  [I/M1/3]  (3.4.1) 

where  In  is  the  natural  logarithm,  I  is  the  impulse  of  the  positive 
portion  of  the  pressure  wave  in  psi-msec  and  M  is  the  mass  of  the 
fish  in  grams.  [The  arbitrary  constant,  80,  was  calculated  by  the 
author  to  correspond  to  the  value,  100,  in  the  Bladder  Oscillation 
Parameter  (Equation  3.1.5).] 

The  results  reported  by  Yelverton  et  al1*  give  a  kill- 
probability 

1 
p  =  (3.4.2) 


1  +  EXP[-0.083(ZI-132) ] 


where  Z   is  the  Impulse  Parameter.   Equation  3.4.2  is  derived  from 
Figure  6  of  Yelverton  et  al  by  fitting  the  cumulative  logistic 
probability  distribution  to  their  results  for  observed  mortality  upon 
holding  the  fish  for  24  hours  after  the  explosion  test. 

Figure  3.4.1a  shows  the  level  3  injury  results  for  Spot  and 
White  Perch  (Tables  3.1.3  and  3.1.4)  plotted  as  a  function  of  the 
Impulse  Damage  Parameter.   Note  that  the  Impulse  Damage  Parameter 
does  not  describe  these  results . 

Figure  3.4.1b  shows  a  subset  of  the  level  3  injury  results 
for  Spot  and  White  Perch  plotted  as  a  function  of  the  Impulse  Damage 


4.   Yelverton,  J.  T.,  et  al,  "The  Relationship  between  Fish  Size  and  their  Response 
to  Underwater  Blast,"  DNA  Report  3677T,  Lovelace  Foundation,  1975. 


47 


NSWC/WOL  TR  76-155 


100% 


(a)  50 


ALL  OF  DATA 


O  OD 

o    o 


o     o  o  o  o     o 


O       O  O  DO 


DO  O  O 


o  o 


IMPULSE  DAMAGE  PARAMETER 


(b) 


100% 


50     - 


IMPULSE  DAMAGE  PARAMETER 


100% 


(c) 


250 


BLADDER  OSCILLATION  PARAMETER,   Z 

FIG.  3.4.1   CORRELATIONS  OF  LEVEL  3  INJURIES  FOR  SPOT  AND  WHITE  PERCH  WITH  IMPULSE 
DAMAGE  PARAMETER  AND  WITH  BLADDER  OSCILLATION  PARAMETER 


48 


NSWC/WOL  TR  7  6-15  5 

Parameter.   In  Figure  3.4.1b  those  data  points  where  the  bladder 
response  calculation  showed  more  than  one  oscillation  during  the  posi- 
tive phase  (Tables  3.1.3  and  3.1.4)  have  been  removed.   The  curve  is 
Equation  3.4.2,  the  kill-probability  result  from  Yelverton  et  al 
(1975) .   Thus,  their  24-hr  mortality  is  approximately  equivalent  to 
our  level  3  dissection  injuries;*  and,  perhaps  more  important,  the 
Impulse  Damage  Parameter  is  shown  to  describe  only  those  explosion 


geometries 

where 

the 

charge 

and/or  the  fish 

are  at  shallow 

enough 

depth 

that 
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ximately  one 

or 

less  cycles 
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face 

reflection 
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the 

positive  pressure 

phase. 

Figure  3.4.1c  shows  for  comparison  all  of  the  level  3  injury 
results  plotted  as  a  function  of  the  Bladder  Oscillation  Parameter. 
The  curve  is  Equation  3.2.1  derived  from  the  Spot  level  3  injury 
data. 

3.5   DISCUSSION  OF  INJURY  CORRELATIONS 

It  is  not  accidental  that  both  the  bladder  oscillation 
parameter  and  the  impulse  damage  parameter  can  be  used  to  describe 
the  fish-kill  for  shallow  explosion  geometries.   To  see  how  this 
comes  about  we  examine  the  response  of  an  air  bubble  to  impulsive 
pressure  loading. 

"Let  the  pressure  be  applied  suddenly  and  let  it  disappear 
again  before  the  bubble  has  had  time  to  change  appreciably  in  size. 
Then,  the  bubble  will  begin  contracting  inward  at  a  certain  radial 
velocity  v.  given  by 


v.  = —  (3.5.1) 

pAi 


*In  examining  dead  and  disabled  fish  which  have  been  collected  from  the  water 

surface  following  underwater  explosions  Martin  Wiley  and  Greig  Peters  (Chesapeake 

Biological  Laboratory,  Solomons,  Md.)  have  found  only  fish  having  injuries  of 

level  3  or  greater.   This  indicates  that  fish  having  received  lesser  injuries  do 

not  show  up  in  the  visible  fish  kill  although  they  may  later  fall  victim  to 
predation. 
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where  p  is  the  density  of  the  water  and  I  =  Jp   dt,  the  applied 
impulse."5   By  the  methods  of  Appendix  A  we  can  now  calculate  the 
radius  ratio  AMAX/AMIN  for  the  oscillating  bubble.   See  Figure  3.5.1. 


The  solution  --  by  means  of  a  tabulated  function  --  is  of  the 


form 


AMAX 
AMIN 


=  FUNCTION 


A.^fpp- 


(3.5.2) 


where  p.  is  the  initial  value  of  the  ambient  water  pressure. 

Since  ^1  p .   does  not  vary  greatly  for  shallow  explosion 
geometries  and  since  in  most  bladder  fish  the  swim  bladder  comprises 
a  roughly  constant  fraction  of  the  total  volume  (about  6%) ,  Equation 
3.5.2  can  be  written 


AMAX 
AMIN 


=  FUNCTION 


1/3 


M 


( approximately) 


(3.5.3) 


where  we  have  substituted  M,  the  mass  of  the  fish,  for  the  volume, 
since  all  fish  are  approximately  neutrally  buoyant.   Eauation  3.5.3 
shows  that  the  Bladder  Oscillation  Parameter  and  the  Impulse  Damage 
Parameter  are  for  practical  purposes  equivalent,  for  the  special 
condition  of  shallow  fish  depth  and  impulsive  pressure  loading. 

Taken  together,  the  results  of  the  Lovelace  Foundation  for 
fish-kill  as  a  function  of  the  impulse  and  the  present  results 
described  in  terms  of  the  bladder  oscillation  parameter  (and,  in 
part  by  the  impulse  damage  parameter)  give  us  confidence  that  we 
have  achieved  a  correct  approximate  solution  to  the  problem  of 


5.   Kennard,  E.  H.,  "Radial  Motion  of  Water  Surrounding  a  Sphere  of  Gas  in 
Relation  to  Pressure  Waves,"  1943,  published  in  Vol.  II  of  "Underwater  Explosion 
Research,"  Office  of  Naval  Research,  1950- 
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predicting   fish-kill.*      The   non-ideal   appearance   of   the    injury 
correlations    shown    in   Figure    3.4.1b   and   c   are   no   doubt   due    to   the 
approximate   nature   of   these   two   damage   parameters.      That   the    impulse 
damage   parameter   can  be   used   to   describe   observed   injuries    for 
pressure   waves   of  more   than    1/2-cycle    in   duration    (Figure    3.4.1b)    is 
clear   evidence   of   its    approximate   nature.      Apparently,    for   shallow 
explosion   geometries    these   two   damage   parameters    approximate    some 
unknown    "true"    damage    parameter   about   equally  well. 

Additional   work   could   possibly  yield   a   better   damage   parameter 
than   either   the   Impulse   Damage   Parameter   or   the   Bladder   Oscillation 
Parameter.      However   such   additional   work   does   not   seem   justified 
because   at   present   our   ability   to   predict    fish-kill    is    limited  by 
other  more   critical    factors,    namely,    knowledge   of   the   pressure 
signature   in   the   presence   of   surface   and  bottom  effects    and   knowledge 
of    the    fish   species/size/density   distribution   in   the   vicinity   of   the 
explosion. 

Finally,    two  observations   are   noted: 

(1)  The   bladder   oscillation    solution   can   be   used   to   define 
the   useful   range  of   the    impulse   damage   parameter. 

(2)  For   shallow  explosion   geometries    fish-kill   computations 
using   the    simpler   impulse   damage   parameter   can   be   used   to   check 
computations   using   the   computed  bladder   oscillation   parameter. 


*We  are  also   confident    (based   on  unpublished  work   to   date)    that    the  Lovelace 
Foundation  results    (Yelverton  et  al,    1975)    can  also  be   described   in   terms   of 
the  bladder  oscillation  parameter.      To   do   so   is   not   a   routine  matter,    however, 
since   the  negative   phase  pressure  and  duration  at   each   fish   location  must  be 
calculated,    and  an  effective  bladder   radius    coefficient,    (Ai)0/L    (analogous    to 
Equations    3.1.7   and    3.1.8),    must  be   determined  by   trial-and-error   calculations 
for  each   fish. 
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4         APPLICATIONS 

4 . 1    FISH-KILL    CONTOURS 

Figure    4.1.1    shows    the   region    of   greater    than    50%    kill    pre- 
dicted   for    18-cm    long    Spot    in    the   vicinity   of   a    32    kilogram   pentolite 
explosion    at    9   meter   depth.      The    procedure    for   making   Figure    4.1.1    is 
as    follows.       First,    we    calculate    the    bladder   oscillation   parameter    Z 
for   an   array   of    horizontal    ranges,    x,    and   depths,    y;    i.e.,    Z    =    f(x,y). 
These    values,    Z    =    100    In (AMAX/AMIN) ,    are    tabulated    in   Table    4.1.1. 
Next,    assuming   that   injury   of   level    3   constitutes   kill,      we    read    from 
Figure    3.2.2    the   value,    Z,-no    =    120,    for   50%   kill   probability.    Finally, 
we    plot    in    Figure   4.1.1    the    regions    of    Table    4.1.1    having    Z    areater 
than    Z50%. 

To   calculate    Z    in   Table    4.1.1   we    first   calculated   an   approxi- 
mate  pressure-time    signature   by   the   method   described    in  Appendix  C. 
Then,    we   calculated   the    response    to   this    signature   as   described   in 
Section    2    and    in  Appendix   B,    except   that   these   computations   were    for 
undamped   oscillatory  bladder   response.* 

Figure    4.1.2    shows   predicted   regions   of   greater   than   10%,    50% 
and   90%    kill    for    21.5-cm   long  White    Perch    for    32    KG   pentolite    at    9 
meter   depth.      The   procedure    for   making   Figure    4.1.2   was    the    same   as 
described    for   Figure    4.1.1.      Values   of    Z,-.0    =    77   and   Zorio    =    190    read 
from  Figure    3.2.2   were   used   to   obtain   the   10%    and   90%   kill    contours 
in    Figure    4.1.2.      In    Figure    4.1.2    contour   details    such    as    shown    in 
Figure    4.1.1   have   been   smoothed  out. 

Note   two   important    features   of   our   solution   to   the    fish-kill 
problem  which   are    evident    in    Figure    4.1.2: 

(1)    The    fish-kill    is    strongly   dependent   on   the   depth   of 
the    fish. 


*Due    to    "no   damping"   and  because   of   discrepancies  between   the  negative   pressures 
and   durations    calculated  by  Appendix  C   and   results    from   the   1975    test   series 
these   calculations   are  considered  approximate   —  but    quite   adequate   for   the 
purposes   of   this   report. 
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TABLE    4.1.1 


BLADDER  OSCILLATION  PARAMETER  Z  AS  A  FUNCTION  OF  FISHES '  DEPTH 

AND  HORIZONTAL  RANGE* 


Charge:   32  KG  Pentolite  at  9  M  Depth 
Fish:   Spot,  18  cm  Fork  Length 


61 


HORIZONTAL  RANGE  (METERS) 
91     122       152      183 


0 

- 

- 

- 

1.5 

29  9 

333 

131 

3.0 

302 

291 

139 

4.6 

310 

288 

111 

6.1 

228 

217 

194 

7.6 

244 

185 

193 

9.1 

174 

207 

150 

10.7 

228 

151 

135 

12.2 

156 

140 

149 

13.7 

205 

159 

140 

15.2 

149 

116 

95 

16.8 

188 

131 

102 

18.3 

142 

138 

132 

19.8 

175 

107 

91 

21.3 

136 

140 

86 

22.9 

163 

95 

106 

24.4 

133 

91 

92 

25.9 

110 

114 

80 

27.4 

128 

96 

95 

29.0 

125 

119 

82 

30.5 

118 

92 

74 

213 


244 


94 

52 

44 

39 

153 

111 

112 

72 

96 

131 

143 

98 

139 

102 

98 

122 

158 

- 

124 

128 

171 

131 

104 

119 

147 

137 

103 

98 

115 

136 

111 

89 

94 

109 

111 

90 

110 

92 

112 

99 

123 

75 

90 

100 

98 

90 

66 

87 

74 

103 

62 

67 

78 

90 

76 

55 

104 

68 

87 

60 

82 

59 

76 

65 

67 

71 

58 

75 

76 

85 

50 

66 

91 

57 

53 

50 

60 

54 

73 

42 

*These  values  for  Z  were  calculated  using  approximate  predicted  pressure-time 
inputs  and  no  damping  of  the  calculated  bladder  motion. 
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(2)  Regions  of  less  than  50%  kill  —  since  they  encompass 
a  tremendous  volume  of  water  —  make  a  major  contribution  to  the 
total  number  of  fish  killed  by  an  underwater  explosion. 

These  two  features  are  examined  further  in  Sections  4.2  and  4.3. 

4.2  KILL  PROBABILITES  AS  A  FUNCTION  OF  DEPTH  FOR  DIFFERENT 
SIZED  FISH 

Figure  4.2.1  shows  the  variation  of  predicted  kill-probability 
with  fishes'  depth  at  a  fixed  horizontal  range  for  different  fish 
corresponding  to  three  different  sized  swim  bladders.   The  predicted 
kill-probability  is  obtained  from  calculated  Z-values,  e.g., 
Table  4.1.1,  using  Equation  3.2.1  (dashed  line  in  Figure  3.2.2).   The 
approximate  length  of  Striped  Bass  corresponding  to  the  equivalent 
bladder  radii  of  the  bladder  oscillation  calculations  was  calculated 
from 

Bladder  Radius  =  >Q42  (4.2.1)* 

Fish  Length 

For  the  two  larger  sized  fish  in  Figure  4.2.1  the  maxima 
occurring  at  about  7  and  11  meters  depth,  respectively,  are  caused  by 
surface  cut-off  occurring  at  the  first  half-cycle  of  bladder 
oscillation.   At  shallower  depths  the  bladder  does  not  have  time  to 
respond  fully  to  the  positive  portion  of  the  explosion  wave;  while 
coincidence  of  the  first  half-cycle  and  surface  cut-off  amounts  to 
a  resonance  between  the  oscillatory  response  and  the  driving  pressure 
field  outside.   Thus,  at  shallow  depth  the  larger  fish  are   in  effect 
protected  from  harm  by  their  swim  bladders;  while  at  the  resonance 
depth  their  swim  bladders  "do  them  in." 


*The  value  .042  is  80%  of  that  experimentally  measured  from  a  sample  of  3 
Stripped  Bass.   The  value  80%  represents  a  guess  at  the  correction  to  go 
from  measured  radius  to  equivalent  computation  radius  for  this  fish. 
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4.3   TOTAL  FISH-KILL  AND  ITS  SPATIAL  DISTRIBUTION 

In  this  section  we  calculate  the  total  fish-kill  from  an 
underwater  explosion  and  its  distribution  as  a  function  of  horizontal 
range  and  depth,  assuming  a  uniform  spatial  distribution  of  a  single 
size  class  of  one  species  of  fish. 

Table  4.3.1  (top)  lists  the  bladder  oscillation  parameter  Z 
for  21.5  cm  White  Perch  calculated  over  a  cylindrical  grid  surround- 
ing a  32  kilogram  charge  fired  at  9  meter  depth.   Table  4.3.1 
(bottom)  lists  the  corresponding  kill  probabilities  from  Equation 
3.2.1. 

To  calculate  the  fish-kill  we  must  know  which  fish  are 
where;  and,  in  general,  this  is  impossible.   Nevertheless,  to  gain 
insight  as  to  the  potential  of  a  given  explosion  geometry  for  killing 
fish  it  is  useful  to  assume  a  fictitious  fish  density  distribution. 
For  this  purpose  we  will  assume  a  nominal  uniform  density  distribution 
of  1  fish  per  thousand  cubic  meters  of  water. 

To  calculate  the  fish-kill  we  must  integrate  the  product, 
fish  density  x  kill  probability,  over  the  region  surrounding  the 
explosion.   Thus,  for  this  example  we  compute  the  products  of  "fish 
present"  and  "kill  probability"  (Tables  4.3.2  and  4.3.1)  to  get  the 
"fish  killed"  in  each  mesh  volume,  Table  4.3.2  (bottom)  —  each  mesh 
consisting  of  a  cylindrical  annulus  10  0  feet  thick  and  5  feet  in  depth. 
Summing  the  distribution  of  "fish  killed"  shown  at  the  bottom  of 
Table  4.3.2  results  in  a  total  kill  of  1245  fish  out  of  6228  fish 
present.   Note  that  we  have  neglected  200  fish  present  in  the  150-ft 
radius  cylinder  of  water  containing  the  charge  as  well  as  all  fish 
present  beyond  a  horizontal  range  of  850  feet.   Nevertheless,  we 
have  described  the  most  significant  portion  of  the  distribution  of 
fish  killed  by  this  underwater  explosion. 

It  is  apparent,  for  example  that  if  fish  were  present  only 
at  depths  greater  than  70  or  80  feet  relatively  few  fish  would  be 
killed.   Likewise,  if  fish  (of  this  particular  size)  were  present 
only  within  some  10  feet  of  the  surface  the  kill  would  be  relatively 
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TABLE  4.3.1         BLADDER  OSCILLATION  PARAMETER  AND  KILL  PROBABILITY  AS  A 
FUNCTION  OF  FISHES'  DEPTH  AND  HORIZONTAL  RANGE 

Charge:    32  KG  Pentolite  at  9  M  Depth 
Fish:    White  Perch,  21.5  cm  Fork  Length 


Fish 

Depth 

(feet) 


Horizontal  Range  (feet) 
200      300   400   500   600   700 


800 


5  - 

251 

208 

41 

41 

35 

10 

7 

10  - 

220 

211 

119 

73 

32 

30 

28 

15  - 

206 

215 

156 

100 

62 

31 

26 

20  - 

287 

197 

131 

129 

86 

79 

58 

25  - 

206 

180 

147 

114 

109 

74 

69 

30  - 

168 

175 

130 

146 

94 

93 

64 

35  - 

172 

187 

117 

120 

113 

88 

81 

40  - 

211 

163 

120 

97 

102 

97 

96 

45  - 

168 

125 

135 

89 

94 

99 

84 

50  - 

144 

110 

139 

94 

76 

8? 

87 

55  - 

133 

111 

121 

99 

70 

77 

73 

60  - 

177 

133 

91 

104 

76 

62 

68 

65  - 

139 

141 

83 

108 

72 

58 

55 

70  - 

132 

110 

81 

94 

86 

63 

51 

75  - 

128 

98 

79 

70 

91 

60 

49 

»0  - 

154 

86 

98 

67 

87 

66 

54 

85  - 

107 

110 

105 

63 

76 

70 

52 

90  - 

120 

119 

81 

61 

56 

75 

57 

95  - 

125 

76 

66 

69 

54 

65 

61 

100  - 

122 

88 

71 

84 

51 

56 

65 

i  BLADDER 

OSCILLATION 
PARAMETER 


200   300   400   500   600   700   800 


Fish 

Depth 

(feet) 


5 

_ 

1.00 

.99 

.01 

.01 

.01 

.00 

.00 

10 

- 

.99 

.99 

.42 

.05 

.01 

.01 

.00 

15 

- 

.99 

.99 

.85 

.20 

.03 

.01 

.00 

20 

- 

1.00 

.98 

.58 

.55 

.10 

.07 

.02 

25 

- 

.99 

.95 

.77 

.35 

.29 

.06 

.04 

30 

- 

.91 

.94 

.57 

.76 

.15 

.15 

.03 

35 

- 

.93 

.97 

.39 

.43 

.34 

.12 

.08 

40 

- 

.99 

.89 

.43 

.18 

.22 

.18 

.17 

45 

- 

.91 

.50 

.63 

.12 

.15 

.19 

.09 

50 

- 

.74 

.30 

.68 

.15 

.06 

.09 

.11 

55 

- 

.61 

.32 

.45 

.19 

.05 

.07 

.05 

60 

- 

.95 

.61 

.13 

.24 

.06 

.03 

.04 

65 

- 

.68 

.71 

.09 

.28 

.05 

.02 

.02 

70 

- 

.60 

.30 

.08 

.15 

.10 

.03 

.02 

75 

- 

.54 

.18 

.07 

.05 

.13 

.03 

.02 

80 

- 

.83 

.10 

.18 

.04 

.11 

.04 

.02 

85 

- 

.27 

.30 

.25 

.03 

.06 

.05 

.02 

90 

- 

.43 

.42 

.06 

.03 

.02 

.06 

.02 

95 

- 

.50 

.06 

.04 

.04 

.02 

.04 

.03 

100 

- 

.46 

.12 

.05 

.09 

.02 

.02 

.04 

KILL 
PROBABILITY 
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TABLE  4.3.2         FISH  PRESENT  AND  NUMBER  KILLED  AS  A  FUNCTION  OF  WATER  DEPTH 
AND  HORIZONTAL  RANGE 

Charge:       32    KG   Pentolite   at    9   M  Depth 

Fish:      White   Perch,    21.5   cm  Fork  Length 


Fish  Density 
Distribution: 


_o  3 

Uniform  Density   =   1   x   10        Fish/Meter 


Fish 

Depth 

(feet) 


Horizontal  Range  (feet) 
200   300   400   500   600   700   800 


5  - 

18 

27 

36 

44 

53 

62 

71 

10  - 

18 

27 

36 

44 

53 

62 

71 

15  - 

18 

27 

36 

44 

53 

62 

71 

20  - 

18 

27 

36 

44 

53 

62 

71 

25  - 

18 

27 

36 

44 

53 

62 

71 

30  - 

18 

27 

36 

44 

53 

62 

71 

35  - 

18 

27 

36 

44 

53 

62 

71 

40  - 

18 

27 

36 

44 

53 

62 

71 

45  - 

18 

27 

36 

44 

53 

62 

71 

50  - 

18 

27 

36 

44 

53 

62 

71 

55  - 

18 

27 

36 

44 

53 

62 

71 

60  - 

18 

27 

36 

44 

53 

62 

71 

6S  - 

18 

?7 

36 

44 

53 

62 

71 

70  - 

18 

27 

36 

44 

53 

62 

71 

75  - 

18 

27 

36 

44 

53 

62 

71 

80  - 

18 

27 

36 

44 

53 

62 

71 

85  - 

18 

27 

36 

44 

53 

62 

71 

90  - 

18 

27 

36 

44 

53 

62 

71 

95  - 

18 

27 

36 

44 

53 

62 

71 

100  - 

18 

27 

36 

44 

53 

62 

71 

FISH 
PRESENT 


200   300   400   500   600   700   800 


Fish 

Depth 

(feet) 


5  - 

18 

26 

0 

0 

0 

0 

0 

10  - 

18 

26 

15 

2 

0 

0 

0 

15  - 

18 

27 

30 

9 

2 

0 

0 

20  - 

18 

26 

21 

25 

6 

5 

2 

25  - 

18 

25 

27 

16 

16 

4 

3 

30  - 

16 

25 

20 

34 

8 

9 

2 

35  - 

17 

26 

14 

19 

18 

7 

6 

40  - 

18 

24 

15 

8 

12 

11 

12 

45  - 

16 

13 

23 

5 

8 

12 

7 

50  - 

13 

8 

24 

7 

3 

5 

8 

55  - 

11 

8 

16 

9 

2 

4 

4 

60  - 

17 

16 

5 

11 

3 

2 

3 

65  - 

12 

19 

3 

13 

3 

2 

70  - 

11 

8 

3 

7 

6 

2 

75  - 

10 

5 

3 

2 

7 

2 

80  - 

15 

3 

7 

2 

6 

2 

85  - 

5 

8 

9 

1 

3 

3 

90  - 

8 

11 

3 

1 

1 

4 

2 

95  - 

9 

2 

1 

2 

1 

2 

2 

100  - 

8 

3 

2 

4 

1 

1 

3 

FISH 
KILLED 
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light.   However,  for  fish  of  this  size  between  depths  of  20  and  60 
feet  the  kill  will  be  extensive  out  to  a  range  of  about  800  feet; 
and  if  a  large  school  several  hundred  feet  in  diameter  with  a  density 
say,  of  one  fish  per  cubic  meter  were  within  this  region,  a 
disastrous  fish-kill  would  result.   (Such^ an- occurrence  is  unlikely, 
however,  as  it  is  Navy  policy  to  delay  operations  if  a  school  of 
fish  is  observed  in  the  vicinity  of  the  test  site.) 
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APPENDIX  A 

RESPONSE  OF  GAS  BLADDER  TO  UNDERWATER  EXPLOSION  PRESSURES: 

SQUARE  STEP  APPROXIMATION6 

ABSTRACT 

This  appendix  presents  the  theory  of  the  radial  motion  of 
a  gas-filled  spherical  bubble  responding  to  step  changes  in  the  water 
pressure. 

The  theory  is  used  to  estimate  the  motion  of  a  fish-bladder 
sized  bubble  subjected  to  the  pressure  signature  from  an  underwater 
explosion.   Some  important  qualitative  results  as  well  as  some  crude 
quantitative  results  are  readily  obtained.   For  example: 

(1)  For  a  simple  positive  step  which  returns  to  ambient 
pressure,  the  maximum  bladder  extension  depends  on  both  the 
magnitude  and  duration  of  the  pressure  pulse.   If  cut-off  occurs 
when  the  bladder  (undamped  oscillation)  is  at  maximum  size,  no 
extension  beyond  its  initial  at-rest  size  occurs.   If  cut-off  occurs 
at  any  other  time,  some  over-extension  of  the  bladder  occurs. 

(2)  A  sudden  decrease  in  the  ambient  pressure  will  always 
result  in  some  over-extension  of  the  bladder.   If  this  neqative 
pressure  (below  ambient)  follows  a  positive  step  with  cut-off 
occurring  at  maximum  size,  the  over-extension  is  due  solely  to  the 
negative  portion  of  the  pressure  signature;  if  not,  the  over-extension 
is  due  to  both  the  positive  step  and  the  negative  pressure  following 
cut-off.   The  effects  are  additive  in  a  non-linear  manner  described 

in  this  appendix. 


6.   Appendix  A  was  formerly  given  internal  distribution  as  NOLTN  10205  ,  Goertner  , 
1974,  "Response  of  Air  Bubbles  to  Underwater  Explosion  Pressures:   Square  Step 
Approximation  with  Application  to  Fish-bladder -cavity  Response". 
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I.   ON  THE  MOTION  OF  AN  AIR  BUBBLE  SUBJECTED  TO  A  SERIES  OF  STEP 
CHANGES  IN  THE  OUTSIDE  WATER  PRESSURE 


1.  Consider  a  spherical  bubble  at  rest  in  water  at  the 
initial  value  p.  of  the  ambient  pressure.   Let  A±   be  its  initial 
radius  when  at  rest.   Its  internal  air  pressure  is  the  initial 
ambient  pressure  p^ .   At  t  =  o,  the  ambient  pressure  jumps  to  a  new 
value  p  ;  and  the  bubble  is  then  oscillating.   The  situation  for  the 
air  pressure  inside  the  bubble,  trying  to  adjust  to  the  new  outside 
pressure,  is  sketched  in  Figure  A-l. 

2.  In  this  Appendix  we  concern  ourselves  with  the  motion 
of  a  bubble  subjected  to  a  succession  of  jumps  to  new  ambient 
pressures  pn.   At  each  ambient  pressure  the  bubble  motion  is 
described  by  the  following  differential  equation  J 


|  (^A3)  A2  +^P0A3  +  E(A)  =  Y  (Al) 


Kinetic  Energy    Potential 
of  Energy  of 

Surrounding  Water   Pushed-back 

Water 


t 

Total    Energy 


Internal   Energy 
of  Air 
Inside   Bubble 


7.      Cole,   R.    H.  ,   19U8,    "Underwater  Explosions",   Princeton  University  Press. 
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m 


M 


Pi  = 

Pl  = 

PM  = 
P  = 


Initial  outside  water  pressure 
New  value  of  outside  water  pressure 
Bubble  air  pressure  at  maximum  size 
Bubble  air  pressure  at  "Equilibrium" 

Pn 


P  =  Bubble  air  pressure  at  minimum  size 


m 


FIG.  A-1      BUBBLE  AIR  PRESSURE  AND  OUTSIDE  WATER  PRESSURE  AS  A 
FUNCTION  OF  TIME  (SKETCH) 
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_,_.   P'(A)  "X  V(A)  ,,„> 

where  E  (A)  =— —+ —  1  (A2) 

A  =  Bubble  Radius 

dA 
A  =  -g^r   ,  the  radial  velocity 

P (A)  =  Air  Pressure  Inside  the  Bubble 
V(A)  =  Bubble  Volume  =  -j^  A 

Y  =  Exponent  in  Adiabatic  P-V  Relationship  =  1.4 

(for  air) 
p  =  Water  Density  =  constant 
P0  =  Ambient  Water  Pressure  =  constant 

Y  =  Total  Energy  of  Bubble  Oscillation  =  constant 

3.   The  oscillatory  motion  described  by  Equation  Al  has 
much  in  common  with  the  oscillation  of  a  simple  mass-spring  system  — 
the  difference  is  that  for  the  bubble  oscillation  the  "mass"  changes 
with  time  and  the  "spring"  is  non-linear.   The  overall  behavior  of 
oscillating  gas  bubbles  was  worked  out  in  considerable  detail  by 
Snay  and  Christian. 

They  rewrote  Equation  Al  in  the  non-dimensional  form 

a3  a2  +  a3  +  k  a"3^"1^!  (A3) 

where  a  is  the  non-dimensional  radius  a=  A/L  and  a  is  the  derivative 

of  "a"  with  respect  to  the  non-dimensional  time  t1  =  t/C .  The 
parameters  k  and  y  characterize  the  bubble  motion.  "k"  can  be 
expressed  by 


k  =    PMm 


Po(Y-D 


i  +  p^M  y 


lm "I 

-i)  J 


where  P„   is  the  internal  pressure  at  either  extremum  of  the  bubble 

Mm  r 

oscillation   —  M  refers    to   the   maximum  bubble   volume    (P      =   minimum 

pressure)  ,    m   refers    to   the   minimum  bubble   volume    (P„   =   maximum 
c  m 

pressure)    —  and   P0    is    the   ambient  hydrostatic   pressure. 


8  Snay,   H.    G. ,   and  E.   A.    Christian,   1952,    "Underwater  Explosion  Phenomena: 
The   Parameters   of  a  Non-Migrating  Bubble   Oscillating  in  an  Incompressible 
Medium",    NAVORD  Report   21*37 
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The  length  scale  factor  L  and  the  time  scale  factor  C 


are  defined  by 


1/2 


(A5) 


(A6) 


where  Y  is  the  total  energy  (Equation  Al)  and  p  is  the  density  of 
the  water. 

5.   Snay  and  Christian  in  their  report  have  determined  the 

overall  behavior  of  the  oscillatory  solutions  to  Equation  A3  for  y 

in  the  range,  1.0  <  y   <1.5.  Table  Al  (below)  list  some  of  their 

results  for  Y  =  1.4.   The  dimensionless  quantities  —  a.,   a   and  T' 

^  Mm 

—  are  the  maximum  radius,  the  minimum  radius  and  the  period  of 
oscillation  of  the  bubble,  respectively. 


TABLE  A-l 
Bubble  Parameters  for  y 


=    1.4 


k 

M 

a  /a,, 
nr  M 

T' 

T  '  Calculated 
From  Equation 

A32 

0.0 

1.000 

0.0 

1.492 

1.476 

0.05 

0.982 

0.085 

1.505 

1.496 

0.10 

0.963 

0.152 

1.522 

1.516 

0.15 

0.942 

0.220 

1.538 

1.536 

0.20 

0.920 

0.288 

1.557 

1.556 

0.25 

0.893 

0.364 

1.577 

1.576 

0.30 

0.863 

0.447 

1.596 

1.596 

0.35 

0.823 

0.548 

1.616 

1.616 

0.40 

0.768 

0.699 

1.636 

1.636 

0.432 

0.659 

1.000 

1.650 

1.649 
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6.  In  what  follows,  we  will  use  both  Equations  Al  and  A3. 

Transitions  from  one  outside  pressure  level,  p  ,  to  another,  p  ,., 

^n  ^n+1 

are  simply  described   using  Equation  Al.  And,  Snay  and  Christian's 
results,  e.g.,  those  listed  in  Table  I,  which  were  computed  from 
Equation  A3,  are  convenient  for  computing  the  bubble  behavior  within 
a  given  state. 

7.  Boundary  Conditions  at  Pressure  Jumps.   We  now  consider 

a  step  change  in  outside  pressure  from  p   to  p  , . .   Let  V   be  the 
c  3  c  ^n    ^n+1        c 

n 

bubble  volume  at  the  instant  of  change,  and  y   and  y  , .  be  the  old 

.,  n     Jn+1 

and  new  total  energies.   For  the  n    state  we  rewrite  Equation  Al  as 


/4tt  p  A3  \a2  +  4iTp   A3  + 

Vs"      )     ^n 


3  /  4tt  p  A   \a   +  4iTp   A   +  E(A)  =Y        (A7) 
2    I  ^J-  — sr  n  n 


Now,  since  the  oscillating  system  described  by  (A7)  has  finite  mass, 
neither  A  nor  A  can  change  impulsively,  i.e., 

AA  =  0  (A8) 

AA  =  0  (A9) 

at  each  pressure  jump.   Thus,  the  first  and  third  terms  of  (A7)  do 
not  change  at  the  jump,  and  the  change  in  the  total  enerqy  is  given 
by 

Yn+1  =  Yn  +  (Pn+l  "  Pn)  Vc  (A10) 

n 


where  V   is  the  bubble  volume  at  the  time  of  the  jump, 
n 

8.   Equations  (A7) ,  (A8),  (A10)  completely  specify  the 

motion  in  the  n    state  given  suitable  initial  conditions.   In  the 

present  case  the  initial  or  zeroth  state  is  specified  by 

A  =  constant  =  A.  (All) 

P(A.)  =  constant  =  p.  (A12) 

l  *i 
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In  particular,  the  total  energy  Y   and  the  outside  pressure  p   are 
sufficient  to  specify  the  amplitude  of  the  oscillation.   And,  if  the 
radius  function  A(t)  is  known,  the  phase  will  then  be  determined  by 
the  jump  conditions,  AA  =  0  and  AA  =  0 . 

+■  v> 
9.   Some  General  Equations  for  the  n    State.   We  now  list 

some  further  equations  that  can  be  used  to  calculate  the  parameters 

+*  v» 
of  the  bubble  oscillation  in  the  n    state. 


kn  =  7T 7r (A13) 

—  (pf-r) 

/3Y    X1/3 

Ln  =  (^wV)  (A14) 


-  I  '  (A15) 


:n  "  Ln(%) 


(AMAXJ 


a  (kn)  x  Ln  (A16) 


AMIX 

AMAX  =    Tabulated   Function    (k    )  (A17) 

n 


'A-  \  3Y 
P(A)    =    p.lp-  J  (A18) 


]/3Y 


(A19) 


T      =    t(k    )     x   C  (A20) 

n  n  n 
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Equation  A19  gives  the  equilibrium  radius  A  (radius  after  the 
oscillation  has  damped  out)  for  the  bubble  in  the  n    state. 
A19  is  derived  from  (A18)  by  setting  P (A)  equal  to  p 


n 

10.   Bubble  Motion  in  State  1.   Consider  the  step  change 
in  outside  pressure  from  p,  to  p..   In  accordance  with  the  jump 
conditions  (A8)  and  (A9)  and  the  initial  conditions  (All)  and  (A12) 
we  set  A  =  0  and  A  =  A.  in  (A7)  to  obtain  the  total  energy 


Yl  =  (Pl  +  Y=T>  Vi 


(A21) 


where  we  have  also  used  Equation  A2  to  evaluate  the  internal  energy 


E(A.) 


11.   Since  A  =  0,  the  oscillation  in  state  1  starts  at  an 


extremum.   Whether  this  extremum  is  a  maximum  volume  V„  or  minimum 

M 

V"m  is  determined  by  p  being  "greater  than"  or  "less  than"  p.. 
If  p,  >  p . ,  then 

AMX  =  Ai  (A22) 


PMi  =  Pi  (A23) 


IF  p.  <  p. ,  then 


\   =   Ai  <A24> 


Pn   -  Pi  (A25) 
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12.   Since  the  oscillation  in  state  1  starts  at  an  extremum, 
the  bubble  parameter  k  can  easily  be  computed  from  Equation  A4 .   From 
(A4)  together  with  (A23)  and  A25  we  get 


kl  =  xi(1+xi)  Y 


(A26) 


where 


xi  ■  Tr^lbr  (A27) 


13.   Calculation  of  the  Subsequent  Extremum.   The  bubble 
radius  at  the  next  extremum  (]/2  period  later  in  time)  is  calculated 
from  the  initial  value  A   =  A.  by  either  dividing  or  multiplying 


by  the  ratio  (  —  J  obtained  from  Table  I  (by  linear  interpolation 

\am/l 
for  the  calculations  done  here) .   The  corresponding  air  pressure  P 

inside  the  bubble  is  then  calculated  using  Equation  A18,  i.e., 


(A28) 


14.   Calculation  of  the  Bubble  Period.   The  bubble  period 
of  oscillation  T  is  the  product  of  the  dimensionless  period, 
t  =  F(y,k) ,  given  in  Table  I  and  the  time  scale  factor  C  (Equations 
A15  and  A14) 


Tx  =  T'  x  C  (A29) 

Combining  (A14,  (A21)  and  (A27) ,  we  write  for  the  length  scale  factor 


Lx  =  Ai  (1  +  X1)]/3  (A30) 
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And,  combining  (A15) ,  (A27)  and  (A30) ,  we  write  the  time  scale 
factor  as 

1/2 
Ci    =    (9^±JL    ^  j       (1  +  Xi)     A_  (A31) 

where  p  is  the  water  density. 

15.   For  our  purposes  the  dimensionless  period  T1,  is 
conveniently  approximated  by  the  linear  function 


T  '  -  1.476  +  0.4  k  (A32> 


Values  of  T  '  calculated  by  Equation  A32  are  listed  in  Table  1 
alongside  the  tabulated  function.   The  approximation,  Equation  A3 2, 
is  good  to  better  than  0.5%  for  values  of  k  }    0.1. 

16.   Finally,  combining  (A26),  (A27) ,  (A29),  (A31)  and  (A32) 
we  write  for  the  period  of  oscillation 


*  "  (t  *) 


1/2  , /. 

(1  +  X]_)  '         (1.476  +  0.4kT)  Ai      (A33) 


The  physical  quantities  in  Equation  A33  may  be  expressed  in  any 
consistent  set  of  units.   For  English  Engineering  Units  we  have 


T,  =  Bubble  Period  (SECONDS) 

3 
p  =  Density  of  the  Water  (SLUGS/FT  ) 

2 
Initial  Value  of  Ambient  Water  Pressure  (LBS/FT  ) 

^i  '   [pressure  in  LBS/FT2  =  (pressure  in  PSI)  x  144] 

A.  =  Initial  Radius  of  Bubble  (FT) 

l 
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17.  Simplifying  Assumption.   To  proceed  further  with  the 
general  case  we  would  now  need  the  bubble  radius  as  a  function  of 
time  in  order  to  apply  the  jump  equations  (A8) ,  (A9)  and  (A10)  at 
the  instant  the  outside  pressure  changes  to  its  new  value,  p- • 
Although  not  particularly  difficult,  such  computations  are  beyond 
the  limited  scope  of  this  study.   Consequently,  at  this  point  we 
introduce  the  simplifying  assumption  that  all  pressure  changes  occur 
at  an  instant  when  A  =  0,  i.e.,  either  at  an  extremum  or  at 
equilibrium  bubble  radius  after  the  oscillation  has  completely 
damped  out.   Pressure  jumps,  p   to  p   , ,  which  happen  to  occur  at 
bubble  extrema  give  the  greatest  and  also  the  least   amplitudes  of 

■f-  V» 

oscillation  in  the  n+1    state  (V   =  V..   in  Equation  A10)  .   Thus, 

c     Mm 

even  when  the  pressure  jumps  occur  at  other  times  than  the  extrema, 
we  can  place  exact  upper  and  lower  limits  on  the  amplitude  of 
oscillation  in  each  state. 

18.  Bubble  Motion  in  State  2.   In  accordance  with  our 

simplifying  assumption  we  will  assume  the  pressure  jump,  p,  to  p_  , 

occurs  at  a  time  t    such  that 

cl 


t    =  N  X  ^  T.  (A34) 

C,  2       1 

where  N  is  any  positive  integer.    Alternatively,  t   may  be  taken 

cl 

large  enough  that  the  bubble  oscillation  has  damped  out,  i.e., 

—  1/3y 

A    =  A,  =  A.  (p,/p.)  '     ',  the  equilibrium  radius. 

19.   Let  V   ,  A   ,  P    be  the  volume,  radius  and  air 
Cl    Cl    °1 
pressure,  respectively,  of  the  bubble  at  time,  t   .   If  the  new  out- 

Cl 
side  pressure  p„  is  higher  than  P   ,  the  air  pressure  inside  the 
l  cl 

bubble,  the  new  oscillation  begins  at  maximum  volume;  and  if  p„  is 

lower  than  P   ,  the  oscillation  begins  at  minimum  volume.   Or, 
c 
1 
restated  we  have 


A-ll 


where 
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If  p  >  P   ,  then    A     =  A  (A35) 

l  C1  M2        C1 


PM    =  P„  (A36) 

M2      C]_ 


If  pn  <  P   ,  then   A     =  A  (A37) 

2     cl  m2      °1 


P     =  P  (A38) 

m2      Cl 


20.   Next,  for  the  bubble  parameter  k  we  have 


k2  =  X2  (1  +  X2)  Y  (A39) 


P 
X.  =  -, 4 (A40) 


L2    (Y-D  P2 


With  this  value  for  k_  the  bubble  radius  at  the  subsequent  extremum 

can  now  be  calculated  by  either  dividing  or  multiplying  the  starting 

value  A   by  the  ratio  (a  AO  »  obtained  from  Table  I.  The  corre- 
c,               m  M  z 

sponding  air  pressure  P    inside  the  bubble  can  then  be  calculated 

by  Equation  A18.   And  finally  the  period  of  oscillation  can  be 

calculated  from 


T„  = 


(§f*2) 


V2 


(1  +X0)]/3   (1.476  +  0.4  k,)  A^       (A41) 


2     ^~   '   2  |      vx  "2'     li,"u     '    2'   c1 
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21.   Bubble  Motion  in  State  n.   The  corresponding  equations 
for  the  n    state  can  easily  be  obtained  from  those  for  state  2 
by  substituting  "n's"  for  '^'s"  and  "n-l's"  for  "l's".   The  only 
apparent  exception  is  that  the  equivalent  equation  for  (A34)  should 
read 


T     =  T      +N  XtT  (A42) 

c       c   .        2  n 
n       n-1 
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II.   ON  THE  MOTION  OF  A  SWIM-BLADDER-SIZE  AIR  BUBBLE  SUBJECTED 
TO  AN  EXPLOSION  GENERATED  PRESSURE  WAVE 


22.  At  the  present  time  we  ask  such  questions  as:   "Can 
a  nearby  explosion  kill  or  seriously  damage  fish  by  extending  their 
swim  bladders?"   Prodded  by  this,  let  us  calculate  the  behavior  of 

a  bladder-sized  air  bubble  subjected  to  an  explosion  pressure 
signature — one  which  kills  many,  but  not  all,  of  the  fish  present. 
The  pressure  signature  shown  in  Figure  A-2  is  derived  from  a  experi- 
mental test  condition  (Shot  525,  Chesapeake  Bay  Tests,  July-Aug  73) 
and  appears  to  meet  this  criterion.   It  corresponds  to  a  fish  cage 
location  where,  upon  subsequent  dissection,  all  ten  of  ten  5.4"  long 
Spot  and  six  of  ten  7.2"  long  White  Perch  were  judged  to  have 
received  lethal  damage.   For  this  calculation  we  will  further 
approximate  the  outside  pressure  signature  with  the  square-stepped 
one  which  is  also  shown  in  Figure  A-2. 

23.  To  get  a  better  feel  for  physical  phenomena,  we  will 
do  the  calculation  in  three  stages.   First,  we  will  calculate  the 
response  to  the  positive  portion,  a  square-step  which  returns  to  the 
initial  ambient  pressure  level;  then,  the  response  to  a  simple  step 
decrease  in  outside  pressure  to  the  level  of  the  explosion  induced 
underpressure;  and  finally,  the  combined  response  to  the  square- 
stepped  signature  sketched  in  Figure  A-2. 

24.  We  will  take  our  bladder-sized  air  bubble  from  a  5.4" 
long  Spot,  and  start  the  calculation  with  an  initial  (at  rest)  bubble 
radius  given  by* 

Aj_  =  0.0395  L  (A43) 


it 
Footnote  on  page  A-16 
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where  L  is  the  fork  length  of  the  fish.   For  our  5.4"  Spot  this  gives 
an  initial  bubble  radius 


A.  =  0.0178  ft  (A44) 


The  initial  ambient  pressure  p^^  =  32  psi,  which  corresponds  to  a  depth 
of  40  feet  in  fresh  water.   The  magnitude  of  the  pressure  step, 
Pl  "  Pi'  is  50  Psi;  and  its  time  duration  TPOS  =  2.3  milliseconds. 
The  period  ^  at  p   =  82  psi  turns  out  to  be  0.569  milliseconds,  so 
that  the  oscillation  in  State  1,  which  starts  at  maximum  bubble  size 


A 

the   time    cut   off   occurs 


M     =  Ai,    goes    through    2.3/0.569    =4.04   periods   of   oscillation   by 


(Footnote    from     page   A-14) 

*This    gives   a    rough   approximation   to   the  bladder   volume   of   a    5.4"    lone 

fish   at    40ft-depth   and    500    feet   distant    from  a    68-lb   pentolite   charge 

(Shot    525) .       Equation   A43   was    conjured   up   by    the    author    for   a 

"nominal    fish".       It   represents    the   radius   of   a    sphere   of    1/20        of 

g 
the   estimated   volume   of   the    fish    (after  Weston) .         The   volume   of   the 

fish   is   taken   as    an   ellipsoid  whose    semi-major   axes    (A,B,C)    bear   the 

relation,    A   =    0.38    L,    B    =    0.30    A,    C    =    0.25    B,    where    L    is    the   overall 

length   of   the    fish.      The    above   coefficients   are    average   values   taken 

from  crude   measurements   of   sketches   of    fish    shown   by  Lippson. 


9.   Weston,   D.    E. ,    1966,    "Sound  Propagation   in  the  Presence  of  Bladder  Fish", 
published  in  Vol.    II   of   "Underwater  Acoustics",    edited  by  V.    M.    Alters, 
1967,   Plenum  Press 

10.    Lippson,   A.    J.,   1973,    "The   Chesapeake  Bay  in  Maryland  —  an  Atlas   of 
Natural  Resources",   John  Hopkins   University  Press 
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25.  At  cut-off  time,  TPOS,  the  measured  outside  water 
pressure  drops  to  about  13  psi  (Figure  A2).   Table  A2  lists  limiting 
case  results  for  the  final  maximum  radius  (in  State  2)  which  were 
calculated  from  the  above  parameter  values  by  assuming  cut-off  times 

TPOS  occurina  at  A.,,  A  ,  and  also  sufficiently  late  that  the 

Mm  ■*      _ 

oscillation  has  damped-out  (so  that  A  =  constant  =  A) . 

26.  Two  sets  of  values  for  A   /A.  are  given  in  Table  A2 , 

M2   1 
one  for  a  final  outside  pressure  p_  =  p.  =32  psi,  the  other  for 

the  measured  value  of   p„  =  13  psi.   The  results  listed  in  the  top 

row  of  Table  A2  (cut-off  at  bubble  maximum)  are  the  same  as  if  the 

positive  outside  overpressure,  p.  -  p.,  had  never  occurred.   In 

other  words,  the  upper  right  hand  value  of  A   /A.  also  applies  to  a 

simple  step  decrease  in  the  ambient  pressure  p.  to  a  new  value, 
p_  =  13  psi. 

27.  It  is  apparent  from  these  crude  results  that  both 
positive  and  negative  excursions  of  the  outside  pressure  significantly 
influence  the  final  overshoot;  and  that  time-of-occurrence  of  cut-off 
relative  to  the  phase  of  the  bubble  oscillation  cycle  is  of  critical 
importance  for  positive  excursions.   It  also  looks  like  a  good  bet 
that  underwater  explosions  do  kill  or  cause  serious  damage  to  nearby 
fish  by  extending  their  swim  bladders. 

28.  Figures  A3  and  A4  sketch  the  qualitative  behavior  of 
the  inside  air  pressure  and  the  bubble  radius  as  functions  of  time 
for  the  six  cases  listed  in  Table  A2 . 
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III.    DISCUSSION 


29.   After  the  above  analysis  was  completed,  George  Young 
brought  to  my  attention  a  similar  analysis  done  in  1943  by 
E.  H.  Kennard  (in  connection  with  the  proposed  use  of  screens  of 
bubbles  as  a  protective  device  against  explosion  shock  waves) .    His 
development  and  presentation  of  the  incompressible  bubble  theory 
(pp.  190-191,  Ref.  5)  is  verv  convenient  for  calculating  bubble 
motions  such  as  the  example  used  in  this  note,  i.e.,  response  to 
jumps  in  constant  outside  pressure  states  restricted  to  those  instants 
when  the  radial  velocity,  A  =  0.   I  redid  the  calculations  summarized 
in  Table  A2  using  Kennard1 s  "bubble  theory"  and  the  bubble  radii  and 
periods  of  oscillation  were  in  all  cases  within  1%  of  my  previous 
values  calculated  using  the  equations  presented  in  this  Appendix. 
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IV.   CONCLUSIONS 


30.   I  hope  that  the  analysis  presented  here  (or  some 
further  development  of  it)  will  prove  useful  in  correlating  observed 
fish  injury  and  kill  with  incident  explosion  pressure  signatures. 
The  present  analysis  is  but  a   crude  first  cut,  however;  and  if 
this  approach  to  understanding  fish  damage  should  appear  fruitful, 
considerable  further  development  will  probably  be  desirable  and 
perhaps  essential  to  achieving  a  useful  correlation.   Nevertheless, 
some  things  are  already  clear.   For  example: 

(1)  By  itself,  a  step  increase  in  the  ambient  pressure 
cannot  extend  the  bubble/bladder  beyond  its  original  size,  since  the 
bubble  oscillates  between  its  original  volume  and  some  minimum  volume 
which  depends  on  the  magnitude  of  the  incident  shock  pressure. 
Extensions  beyond  the  initial  at-rest  volume  can  only  occur  after 
return  to  ambient  pressure  or  below. 

(2)  For  a  simple  positive  step  pressure  which  returns  to 
the  ambient  level,  the  maximum  bubble/bladder  extension  depends  on 
both  the  magnitude  and  duration  of  the  pressure  pulse.   If  cut-off 
occurs  at  the  instant  of  maximum  bubble  volume,  the  bubble  is 
returned  to  its  initial  at-rest  state,  having  undergone  no  extension 
beyond  this  initial  at-rest  state.   If  cut-off  occurs  at  any  other 
time,  some  extension  beyond  the  initial  at-rest  volume  occurs.   For 
a  given  incident  pressure  level,  the  greatest  bubble/bladder 
extensions  occur  when  cut-off  is  at  the  instant  of  minimum  bubble 
volume. * 


*  Footnote  next  page 
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(3)    A   step  decrease    in   the   ambient   pressure   always   results 
in    an    extension   of    the   bubble/bladder.       If    such    a    drop    follows    a 
positive   step  where   cut-off   occurs    at  maximum  bubble   volume,    the 
extension   beyond    initial   at-rest   volume    is   due   solely   to   the   negative 
portion  of   the   pressure   signature;    if   not,    the   over-extension   is   due 
to   both   the   positive    step   and   the   underpressure    following   cut-off. 
The   effects   are   additive    in   the   non-linear   manner   described   in   this 
appendix. 


(Footnote    from  preceding   page) 

*For  a  simple  step  pressure  pulse  of  level  p    ,   the   energy  of  the  residual  bubble 
oscillation   is   given  by 

Y'2  =  X2-Y.    =    (Pl   -  p.)    x    (V.    -  7Ci) 

where  V   is  the  bubble  volume  at  the  instant  of  cut-off.   Thus,  the  energy  of  the 

Cl 

residual  oscillation  and,  consequently,  extension  of  the  bubble/bladder  beyond 

initial  at-rest  volume  take  on  greatest  values  for  V   =  V  and  tend  to  zero  as 

c     m 

V        approaches  V    . 
c        rr  M 
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APPENDIX  B 

METHOD  FOR  CALCULATING  GAS  BLADDER  RESPONSE 
TO  EXPLOSION  PRESSURE  WAVE 

Figure  2.3.2  sketched  the  way  we  approximated  the  explosion 
pressure  signature  by  a  sequence  of  pressure  steps  in  order  to 
calculate  the  bladder  response  by  means  of  the  equations  developed 
in  Appendix  A.   Appendix  B  gives  the  details  of  the  procedure. 

Boundary  Condition  at  Pressure  Jumps.   Step  changes  in  the 

outside  water  pressure  from  p   to  p  , ,  occurring  at  half-period 
c  ^n    ^n+1         ^         c 

intervals  change  the  equilibrium  pressure  P0  of  the  oscillating 
bubble  flow  (equal  to  the  outside  pressure  p  ) .   Since  the  oscillatory 
system  described  by  Equation  Al  has  finite  mass,  neither  A  nor  A  can 
change  impulsively,  i.e., 


A  A  =  0  (Bl) 

and 

A  A  =  0  (B2) 

for  all  pressure  jumps.   And,  thus  for  pressure  jumps  occurring  at 
extrema  where  A=0 ,  only  the  equilibrium  pressure  P„  can  change. 

Let  A   be  the  maximum   or  minimum  bubble  radius  at  the 
c 
n 

time  of  the  pressure  jump  p   to  p  ,  ,  .   If  A    is  greater  than  the 
t-        j   f  t-n    ^n+1       c      3 

n 

new  equilibrium  radius*  A  ,  ,  ,  the  new  oscillation  begins  at  maximum 

n+1 


*  _ 

The  equilibrium  radius  A  is  the  at-rest  radius  which  corresponds  to  the 

equilibrium  pressure  E,  (=  p  ,  the  outside  pressure).   A   is  calculated 

^      n  n 

from  p  using  Equation  B6,  below. 
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size;  and  if  A    is  smaller  than  A  ,.  the  new  oscillation  begins  at 
c  n+1 

n 

minimum  size.   Or,  restated  we  have 


If         A    >  A  ..  ,  then  AMAX  .-  =  A  (B3) 

c     n+1  »  n+1    c 

n  n 


If         A    <  A  ..  .  then  AMIN  ,.    .  (B4) 

c      n+1  >  n+1  =  A 

n  c 

n 


Equilibrium  Bubble  Radius.   At  any  time  during  the  motion 
the  internal  gas  pressure  p  is  related  to  the  bubble  radius  by 


P(A)  =  pi^l)  (B5) 


Setting   P(A)  equal  to  p   in  (B5)  yields  the  equilibrium  bubble 
radius , 

-/3Y 


/Pi\  X 

A   =  A.  I  — 


(B6) 


where  y   is  the  adiabatic  exponent  for  air  (=  1.40)  and  the  subscript 
"i"  refers  to  the  initial  ambient  "at  rest"  state. 

Equilibrium  Bubble  Period.   The  equilibrium  bubble  period 

T   (for  vanishingly  small  oscillations  about  the  equilibrium  radius) 

n    .  5 

was  given  by  Kennard. 


T   =  2ttA   /p/3Vp  (B7) 

n      n   '   rn 


3 
where  p  is  the  water  density  =  1.940  slugs/ft  . 
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Inserting  the  values  for  the  water  density  and  adiabatic  exponent, 
equation  B7  becomes 


_  A 

T   =  29.7— =S—   milliseconds  (B7a) 

n        /^ 


where  A   is  in  inches  and  p   is  pounds  per  square  inch. 

Maximum  and  Minimum  Radius,  and  Period  of  Oscillation. 
With  increasing  amplitude  the  period  increases;  the  ratio  of  the 
periods,  T/T,  is  plotted  as  a  function  of  the  maximum  radius  AMAX 
in  Figure  Bl.   The  same  figure  also  shows  the  minimum  radius  AMIN 
plotted  as  a  function  AMAX.   Both  curves  were  calculated  by  the 
method  given  in  Appendix  A. 

Some  Details  of  the  Computation.   Consider  the  computation 
shown  in  Figure  2.3.2.   We  will  assume  that  the  pressure  signature 
begins  with  a  shock  front  followed  by  an  approximately  exponential 
decay.   To  get  the  average  outside  pressure  for  the  first  half-period 
of  the  motion  we  estimate  the  period  of  oscillation  by 


T  =cxx  T*  (B8) 


where  T*  is  the  equilibrium  period  calculated  using  the  pressure 
value  at  the  start  of  the  motion  (i.e.,  the  Shockwave  peak  pressure  + 
initial  ambient)  and  oc  is  an  arbitrary  constant.   (For  the  present 
computations  an   ©(-value  of  1.25  was  used  for  the  first  half-cycle 
of  oscillation.) 

The  parameter  T*  is  calculated  from  Equation  B7  using  a  corresponding 
A*  calculated  from  Equation  B6 .   The  starting  point  is  the  initial 
bubble  or  bladder  radius  A.  which  was  estimated  from  experimentally 
measured  values.   For  the  present  calculations  the  ratio  (A. ) 0/L 
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(Equations  3.1.7  and  3.1.8)  was  then  varied  to  obtain  the  best 
possible  correlation  between  the  calculated  and  experimentally 
observed  in  juries  —  Section  3.2  and  3.3. 

Next,  we  read  or  calculate  the  outside  pressure  at  the  end 
of  the  first  half-cycle  of  motion  using  the  period  estimate  calculated 
by  (B8)  ;  and  then  calculate  p..  ,  the  constant  outside  pressure  for 
this  approximate  solution,  as  the  average  of  the  pressures  at  the 
beginning  and  end  of  the  half -cycle.  Using  p  we  can  then  calculate 
the  parameters  for  the  first  half-cycle  of  motion;  and,  in  particular, 
obtain  a  considerably  improved  value  for  T.   Using  the  improved 
value  for  T  we  then  calculate  from  Equation  B8  a  better  value  for  c< 
which  we  use  to  calculate  the  next  half-cycle  of  motion.* 

Surface   Cut-off.   Upon  arrival  of  the  reflected  wave  from 
the  water  surface  the  outside  water  pressure  suddenly  drops  to 
below  the  initial  ambient  level.   This  phenomenon  is  known  as 
"surface  cut-off".   When  cut-off  occurred  near  an  extremum — i.e., 
within  0.10  cycles  of  a  maximum  size,  or  within  0.08  cycles  of  a 
minimum — we  calculated  the  final  motion  at  the  lowered  outside 
pressure  as  starting  from  rest  from  the  nearby  extremum.   Interme- 
diate cases  were  calculated  as  starting  from  rest  using  an  interme- 
diate initial  bubble  size — based  on  the  parameter  (At/T)    which 
was  calculated  as  follows: 


If   A      +  AMAX  ,  then  (At/T)-,,  =  (TPOS  -  t     ) /T  (B9) 

c   ,        n  PC  c   ,    n 

n-1  n-1 


*If  necessary,  we  can  use  this  improved  value  for  T  to  get  a  better  value  of  the 
average  pressure  p  .   So  far  this  has  not  been  necessary. 
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If  A      =  AMIN  ,  then  (At/T)    =  (TPOS  -  t     ) /T   +  1/2     (BIO] 
Cn-1        n  FC  cn-l    n 


The  parameter  (At/T)    locates  positive  phase  cut-off  in 
terms  of  bubble  oscillation  cycles.   "Zero"  is  taken  as  the  last 
expansion  occurring  during  the  positive  phase.   Using  (At/T)  _,  }the 
starting  radius  for  the  final  oscillation  at  lowered  pressure  was 
then  selected  from  the  approximate  values  listed  in  Table  B-l. 

Phase  Shift  Correction.   The  effect  of  the  negative 

pressure  phase  is  often  limited  by  its  short  duration  relative  to 

the  bladder  period  of  oscillation.   Thus,  it  is  also  necessary  to 

account  for  the  phase  shift  which  occurs  upon  surface  cut-off.   This 

was  done  by  adding  a  time  correction  to  the  duration,  At     of  the 
1  ^  neg, 

negative  pressure  phase.   Using  the  parameter  (At/T)    computed 
from  Equation  11  or  12,  this  correction  was  approximated  as  follows: 


If   0  <_(At/T)     £  0.1,  then  Phase  Corr  =  0  (Bll) 


If  0.1  <  (At/T)__  <  0.5, then  Phase  Corr  =  -(t   -  TPOS)     (B12) 
—         PC  —     '  c 

n 


If  0.5  <_    (At/T)    <_   1.0;then  Phase  Corr  =  TPOS  -  t         (B13) 

cn-l 


where  t    is  the  time  at  which  AMIN  would  have  occurred  had  not 

cn 
surface  cut-off  intervened.   (Since  this  subsequent  minimum  depends 
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TABLE  B-l 
APPROXIMATE  BOUNDARY  CONDITION   FOR  USE  WHEN  PRESSURE  JUMPS 

OCCUR  BETWEEN  EXTREMA 


If  (-L 


START   NEXT 

HALF-CYCLE 

0.10 

WITH: 

< 

AMAX 

< 

0.25 

(AMAX   +    A)/2 

< 

0.32 

A 

< 

0.42 

(A    +    AMIN)/2 

< 

0.58 

AMIN 

< 

0.68 

(AMIN    +    A) /2 

< 

0.75 

A 

< 

0.90 

(A    +    AMAX)/2 

< 

1.00 

AMAX 
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primarily  on  the  kinetic  energy  at  cut-off  and  not  on  the  outside 
pressure,  this  t  -value  calculated  using  the  higher  pressure  should 
still  be  a  good  approximation  to  the  time  of  occurrence  of  this 
minimum. ) 

Final  Oscillation  at  Negative  Pressure.   In  calculating 
the  final  bubble  expansion  achieved  during  the  negative  pressure 
phase  we  only  consider  those  cases  in  which  the  final  oscillation 
starts  with  a  minimum  radius--only  they  are  of  interest.   In  order 
to  estimate  this  final  expansion  we  calculate  the  parameter 


(At/T)    =  (DTNEG  +  Phase  Corr)  /  T    +  1/2  (B14) 


where  DTNEG  is  the  duration  of  the  negative  phase — and  again  refer 
to  Table  B-l. 

The  parameter  (At/T)    locates  the  end  of  the  negative 

phase  in  terms  of  bubble  oscillation  cycles.   "Zero"  is  taken  as 

the  last  expansion  occurring  during  the  positive  phase. 

If  (At/T)Np  _<  0.58,  we  conclude  the  negative  phase  is  too  short  for 

any  significant  expansion  to  take  place.   Otherwise,  we  calculate 

the  maximum  size  as  indicated  in  the  table.   And,  if  (At/T)    >  0.90, 

the  maximum  bubble  size  is  taken  equal  to  AMAX.   In  any  event,  if 
(At/T)    <  0.90,  we  also  calculate  the  final  oscillation  that  would 
occur  if  surface  cut-off  had  returned  the  outside  pressure  directly 
to  ambient.   If  this  results  in  a  greater  maximum  expansion,  we 
take  this  value  as  the  maximum  size  achieved  during  final  expansion. 

Damping.   Equation  Al  describes  undamped  radial  pulsation 
of  an  ideal  frictionless  fluid.   While  it  seems  unnecessary  to  make 
a  detailed  study  of  energy  dissipation  in  the  fishes'  oscillating 
swim  bladder,  it  does  seem  prudent  to  at  least  approximately  account 
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for  the  fact  that  such  dissipation  must  occur.   In  this  study  we 
did  so  by  withdrawing  a  fraction  3  of  the  energy  of  the  oscillatory 
motion  at  each  extremum. 

The  total  energy  at  an  extremum  can  be  calculated  from 
Equation  Al  by  dropping  the  first  term  and  substituting  from  (A2) 
and  (B5)  to  get 


n 


Hn    y-1  \   V. 


-Y 


*VMm*n 


(B15) 


4       3 

where  (VM  )   is  the  bubble  volume,  *n  (A..  ) 


Rewriting  (B5)  as 


p(V)  =  pi  (V/V±) 


-Y 


(B16) 


and  noting  that  p   =  p(v  )=  P- (V„/V. )  '  and  substituting 

n       n    ci      n'  i 

into  (B15)  we  get 


Y  =  -JL-  Pn  V 

n    y-1   n   n 


(B17) 


for  the  energy  of  the  non-oscillating  equilibrium  bubble. 

Subtracting  y  from  y  we  get  the  bubble  energy  of 
oscillation , 


Yn    Yn    Yn 


(B18) 
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Finally  combining  (B15) ,  (B17)  and  (B18)  and  again  noting  that 
p   =  p.  (V/V.)  ^  we  get  the  dimensionless  equation, 


Y 


Y-l   ,  1 

Y      Y 


V, 


Mm 
V 


-Y 


V, 


Mm 
V 


-1 


(B19) 


for  the  bubble  energy  of  oscillation. 

At  each  extremum  the  bubble  volume  (or  radius)  was  changed 
to  yield  a  new  energy  of  oscillation, 


Y"  =  (1  -3)  Y1 


(B20) 


For  this  study  we  took 


3  =  0.30 


(B21) 


Considering  the  precision  of  our  present  experimental  data,  the 
precise  value  chosen  for  6  does  not  appear  to  be  critical. 
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APPENDIX  C 
APPROXIMATE  METHOD  FOR  CALCULATING  THE  PRESSURE-TIME  SIGNATURE 


Figure  3.1.1  shows  the  idealized  form  of  pressure-time 
signature  used  for  the  bladder  response  computations  of  this  report. 
For  planning  the  1975  series  of  explosion  tests  and  for  the  examples 
presented  in  Section  4  of  this  report  we  needed  to  predict  for 
arbitrary  explosion  geometries  the  variables,  PMAX,  9,  TPOS ,  PNEG, 
DTNEG,  shown  in  Figure  3.1.1.  This  Appendix  tells  how  this  was  done. 


PRESSURE  SIGNATURE  BEFORE  SURFACE  CUT-OFF.   The  initial 
positive  portion  of  the  signature  is  the  direct  wave  from  the  under- 
water explosion.   The  sudden  negative  excursion  is  caused  by  the 
arrival  of  the  reflected  rarefaction  wave  from  the  water  surface. 
This  phenomenon  is  known  as  "surface  cut-off".   We  calculate  the 
time,  TPOS,  of  this  arrival  relative  to  the  direct  arrival  by 


R..-R 
TPOS  =  — (CI) 


where  R  is  the  slant  range  from  the  charge,  R,  is  the  slant  range 
from  the  image  of  the  charge  mirrored  in  the  water  surface,  and  c0 
is  the  sound  speed  in  the  water.   See  Figure  CI. 

We  calculate  the  free-field  portion  of  the  signature  using 
a  modification  of  the  empirical  Shockwave  Similitude  Equations.   The 
peak  pressure  PMAX  and  the  decay  constant  8  are  calculated  by  the 
usual  similitude  relations 
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l/3\a 
PMAX  =  k  (L-  (C2) 


(^7 


a.  w 


V3  /_R    \P  (C3) 


O3) 


where  W  is  the  explosive  weight  of  the  charge,  R  is  the  slant  range 
from  the  charge,  and  k,  a  ,  £,  and  8  are  empirical  constants  (see, 
e.g.,  Reference  7).  Values  for  the  constants,  k,  and  I   and  $  used 
for  the  calculations  of  this  report  can  be  obtained  from  Equations 
3.1.1  and  3.1.3. 

The  instantaneous  free-field  (or  direct  arrival)  pressure, 
p_,  was  calculated  using  two  separate  exponential  segments  joined 
at  t  =  1.8  0 


PMAX  e"t/G  (t  <  1.8  0)  (C4a) 

PD(t)=j 

0.25  PMAX  e"t/4'3  6        (t  >  1.8  0)  (C4b) 


which  were  fitted  to  the  pressure  signatures  recorded  on  the  1973 
explosion  test  series. 

CALCULATION  OF  PNEG,  NEGATIVE  PRESSURE  FOLLOWING  SURFACE 
CUT-OFF.   PNEG  was  determined  by  a  calculation  along  the  surface 
reflected  ray  from  the  charge  to  the  fish.   If  the  surface  reflected 
ray  went  thru  the  region  of  bulk  cavitation,  PNEG  was  computed  from 
the  depths  at  which  the  ray  intersected  the  top  and  bottom  of  the 
cavitation.   Otherwise,  PNEG  was  computed  from  the  linear  super- 
position of  the  direct  and  surface-reflected  pressure  waves. 
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Calculation  of  PNEG  Above/Inside/Under  Region  of  Cavitation 
Let  the  depths,  y   and  y  ,  locate  the  top  and  bottom  of  the  region 
of  cavitation.   Let  p   =p.+PNEG  be  the  absolute  pressure  under  the 

Ado   1 

influence  of  cavitation,  i.e.,  following  arrival  of  the  surface 
reflected  wave.   For  fish  located  above/inside/under  the  region  of 
cavitation,  p,  DC.  is  determined  as  follows: 

Ado 


Above:  pAT3C,  changes  linearly  with  depth  from  PATM  at 

Add 

the  surface  to  PVAP  at  the  top  of  the  cavitation. 
Inside:  p,DC.  =  PVAP 

Ado 

Under:   pABg  =  PVAP  +  pg(y-yc) 

where  PVAP  and  p  are  the  vapor  pressure  and  density  of  water,  and 
g  is  the  acceleration  of  gravity.   Since 


PNEG  =  Pabs  -  Pi  <C5) 

the  corresponding  equations  for  PNEG  are  (see  Figure  C2) 

-[PATM  +  pgy   -  PVAP]^—    0<y<y  (C6a) 

Yt 

PNEG  =)     -[PATM  +  pgy  -  PVAP]     yt<y£y  (C6b) 


-[PATM  +  pgy   -  PVAP]      y>y  (C6c) 
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FIG.  C-2  VARIATION  OF  PNEG  ABOVE/INSIDE/UNDER  REGION  OF  BULK  CAVITATION 
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Calculation  of  Depth,  y  ,  of  Bottom  of  Cavitation.   Given 

Equations  C6a,  b,  c  the  problem  of  determining  PNEG  reduces  to 
determination  of  the  depths  y   and  y   of  the  top  and  bottom  of  the 
cavitation.   To  determine  y  we  plotted,  as  a   function  of  the  charge 

weight  W,  values  of  y   calculated  from  PNEG-values  measured  on  the 

c  1 

1973  Chesapeake  Bay  Test  Program    using  Equation  C6c.   The  equation.* 


y   =  4.4  W0,3  (C7) 

c 

where  W  is  the  explosive  weight  in  pounds  pentolite  and  y   is  in 
feet  was  used  to  approximate  these  measurements. 

Calculation  of  Depth,  y  ,  of  Top  of  Cavitation.   The  top 

(beginning)  of  the  cavitation  region  was  taken  as  the  point  on  the 
reflected  ray  where  the  reflected  shock  first  lowered  the  absolute 
pressure  to  the  vapor  pressure  of  the  water.  To  calculate  y  we  use 
the  linear  superposition  of  the  direct  and  surface-reflected  pressure 
waves , 


PSUM(t)  =  PD(t)  +  PR(t)  (C8) 


The   direct  wave,    p    (t) ,    is    given   by   equations   C2 ,    3,    4.      The 
reflected  wave,    p    (t) ,    can   also  be   calcuated   using   Equations   C2,    3,     4 
provided   one   makes    the   substitutions    "-k"    for    "k"    and    "R, "     (slant 
range    from   image)    for    "R".      To   determine   y     we   then    solve    for   the 
point   on   the   reflected   ray   path  where 


*Equation  CT   gives   a   crude   empirical   fit   to  the  PNEG-values  measured  on  the  1973 
test   program.      It  was   used  out   of  necessity  to  make  kill  probability  calculations 
used  to  design  the  1975   test  program  hut  has   not  been  re-evauluated  using  those 
and  other  available  PNEG-measurements. 
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p.  +  pn(TPOS)  +  p_(TPOS)  =  PVAP  (C9) 

1      D  i\ 


Calculation  of  PNEG  when  Reflected  Ray  is  Beyond  Region  of 
Cavitation.   If  the  total  pressure,  p.  +  p  (TPOS)  +  p  (TPOS) , 

ID  R 

0  3 
upon  arrival  of  the  surface  refection  at  depth,  y   =  4.4  W 

(Equation  C7) ,  was  greater  than  or  equal  to  PVAP,  the  reflected  ray 

was  considered  to  lie  beyond  the  region  of  cavitation.   When  this 

happened  PNEG  was  approximated  by 


PNEG  =  |[pD(TPOS)  +  pR(TPOS)]  (CIO) 


calculated  at  the  location  of  the  fish. 

CALCULATION  OF  DTNEG,  DURATION  OF  NEGATIVE  PRESSURE,  PNEG. 
DTNEG  was  also  determined  by  a  calculation  along  the  surface" 
reflected  ray  from  the  charge  to  the  fish.   If  the  surface-reflected 
ray  went  thru  the  region  of  cavitation  and  the  condition, 


R/W1/3  <  80  (Cll) 


where 

R  =  slant  range  from  charge  to  fish  in  feet 
W  =  explosive  weight  (pentolite)  in  pounds 

was  true,  then  DTNEG  was  calculated  from  the  time-of-f light  of  the 
water  layer  on  top  of  the  cavitated  region.   Otherwise,  DTNEG  was 
computed  from  the  linear  superposition  of  the  direct  and  surface- 
reflected  pressure  waves. 
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Calculation  of  DTNEG  from  Time-of-Flight  of  Water  Layer. 
The  duration,  DTNEG,  of  the  negative  phase  corresponds  to  the 
duration  of  bulk  cavitation  in  the  neighborhood  of  the  fish.   For 
this  study  the  duration  of  cavitation  was  calculated  at  the  point  on 
the  reflected  ray  at  depth  y   (closure  depth)  given  by  Equation  C7 . 
Figure  C3  shows  the  geometry  of  the  problem.   To  get  a  first  approxi- 
mation to  the  duration  of  bulk  cavitation  at  this  point  we  used 
Walker  and  Gordon's  result  for  time  of  flight  of  a  water  layer  of 
thickness,  y  ,  decelerating  (falling  back)  due  to  gravity  and  atmo- 
spheric  pressure 


2  •  PMAX  ♦  6 

TFLIGHT  =  —7 .  .C.  _LCnAmM  (C12) 

pg  (y   +  k)  +  PATM 


where  PMAX   and  9   are  the  Shockwave  peak  pressure  and  decay  constant 

at  the  closure  point  (x  ,y  ) , 

and 


k  = 


e  r 

c   z 


DOB 


and 


c0    =   sound   speed    in  water 

R      =    slant    range    from  charge   to   point   of   reflection   at 

water    surface 
DOB   =   depth    of   burst. 


11/  Walker,   R.    R. ,    and  J.    D.    Gordon,    1966,    "A  Study  of  the  Bulk  Cavitation   Caused 
by  Underwater  Explosions",   David  Taylor  Model  Basin  Report   1896,   pl8 , 
Equation   5.      Note,   however,   that   we  used  an   empirical   equation   (CT)    to   calculate 
the  closure   depth  rather  than  Walker   and  Gordon's   result. 
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To  determine  DTNEG,  we  plotted  the  ratio  of  TFLIGHT  (calculated  by 

C12)  to  DTNEG-values  measured  on  the  1973  Chesapeake  Bay  Program   , 

as  a  function  of  R  /AMAX,  where  R   is  the  slant  range  from  charae  to 

c  c 

cavitation  closure  point  and  AMAX  is  the  maximum  radius  of  the 
explosion  bubble.   This  plot  is  represented  by  the  empirical  equation* 


TFLIGHT 


DTNEG  = 


4.04  -  1.13  ln(R  /AMAX) 


(C13) 


where 


AMAX  =  12 .7 


W 


DOB  +  (33.9/14.7) -PATM. 


V3 


(C14) 


and,  AMAX  and  DOB  are  in  feet,  PATM  is  in  pounds  per  square  inch, 
and  W  is  the  explosive  weight  in  pounds  pentolite. 

Thus,  in  summary,  when  the  reflected  ray  intersected  the 
region  of  cavitation  and  Equation  Cll  was  true,  Equations  C12,  13,  14 
were  used  to  calculate  DTNEG. 

Calculation  of  DTNEG  from  Superposition  of  Direct  and 
Surface-Refelcted  Waves.   If  the  reflected  ray  did  not  intersect  the 
region  of  cavitation  or  Equation  Cll  was  false,  DTNEG  was  calculated 
from  the  linear  superposition  of  the  direct  and  reflected  pressure 
waves  by  means  of  an  empirical  equation  adjusted  to  the  1973 
Chesapeake  Bay  Test  data.   This  empirical  equation  was 


*This  empirical  fit  was  used  to  make  the  kill  probability  calculations  for  design- 
ing the  1975  test  program.   It  has  not  been  re-evaluated  using  the  1975  test 
program  and  other  available  DTNEG-measurements . 
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DTNEG  =  p~yM(PNEG/10)  -  TPOS  (C15) 


where  PqnM  is  the  inverse  of  Equation  C8,  TPOS  is  the  positive 
duration,  and  PNEG  is  calculated  by  the  appropriate  equation,  C6  or 
CIO.   Thus,  the  end  of  the  negative  phase  was  taken  as  the  point 
where  the  negative  pressure  calculated  by  Equation  C8  returned  to 
10%  of  the  value,  PNEG. 
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APPENDIX  D 
LOCATION  OF  FISH  CAGES  BY  SOUND  RANGING 


On  four  of  the  six  shots  of  the  1975  tests,  water  currents 
caused  significant  deviation  of  the  charge  and  the  fish  cages  from 
their  intended  positions.   The  charge  was  suspended  on  one  line  and 
the  cages — with  a  pressure  gage  attached  to  each — where  suspended  on 
another.   In  this  appendix,  we  consider  the  problem  of  determining 
the  deviations  of  the  charge  and  the  string  of  fish  cages  from  their 
intended  locations. 

The  separation,  x0 ,  of  the  two  lines  at  the  surface  was 
assumed  known.   Figure  Dl  shows  the  geometry  of  the  problem.   The 
length  of  the  charge  support  line  and  the  lengths  of  line  between 
gages  were  also  assumed  known.   We  further  assumed  that  the  charge 
and  gages  lie  in  the  the  same  vertical  plane. 

We  solved  this  problem  using  the  measured  arrival  times 
of  the  direct  Shockwave  and  its  reflected  tension  wave  from  the 
water  surface  at  the  respective  gage  locations.   Step  1  was  to  use 
the  relative  arrivals  of  the  direct  shock  to  determine  an  orientation 
for  the  gage  line  giving  a  set  of  constant  differences  between  the 
calculated  arrivals  and  the  measured  ones.   Step  2  was  to  use  these 
gage  positions  to  calculate  the  surface  cut-off  time — the  time 
interval  between  arrival  of  the  direct  shock  and  its  reflection-- 
and  to  compare  these  values  with  the  corresponding  measured  values. 
If  there  existed  a  systematic  discrepancy  in  surface  cut-off  times 
not  accountable  to  variations  in  the  sound  speed  over  the  two  ray 
paths,  we  then  adjusted  the  angle,  8  ,  of  the  charge  support  line 
and  redid  steps  1  and  2.   The  calculations  were  redone  until  there 
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was  no  systematic  discrepancy  in  surface  sut-off  time  not  accountable 
to  variation  in  sound  speed  over  the  two  ray  paths. 

The  computations  of  steps  1  and  2  assumed  that  the  sound 
speed  did  not  vary  with  depth.   R.  S.  Price  of  this  center  has 
developed  a  high-speed  computer  program  which  takes  into  account  the 
sound  speed  variation  with  depth--but  assumes  that  the  gages  are  in 
a  straight  line.   For  the  computations  done  here,  these  two  effects-- 
curvature  of  the  gage  line  and  sound  speed  variation--are  the  same 
order  of  magnitude.   In  these  computations,  errors  of  about  0.4 
milliseconds  are  caused  by  variations  in  sound  speed.   These 
correspond  to  gage  position  errors  of  about  2  feet.   For  our  purposes 
errors  of  this  magnitude  were  acceptable,  however,  significantly 
improved  precision  could  be  obtained  by  accounting  for  both  gage  line 
curvature  and  sound  speed  variation  in  the  calculations . 

As  a  final  check  on  the  computations  we  compared  the 

absolute  time  of  arrival  of  the  direct  Shockwave  measured  on  the 

pressure  gage  records  to  that  calculated  for  the  Shockwave  travelling 

between  our  computed  charge  and  gage  locations.   To  do  this  we 

calculated  the  time,  At     ,  between  the  electrical  firing  pulse 

meas  ^  ^ 

and  direct  Shockwave  arrival  as  follows 


At      =  At     ,  +  VEL.  CORR.  +  FIRING  DELAY  (Dl) 

meas      sound 


where 


At     ,  =  transit  time  for  a  sound  wave 
sound 

VEL.  CORR.  =  correction  to  At     ,  to  account  for  detonation 

sound 

wave  velocity  inside  the  charge  and  Shockwave 
velocity  in  the  water 
FIRING  DELAY  =  dwell  time  between  firing  pulse  and  initiation 
of  charge 
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The  firing  delay  depends  on  the  firing  circuit,  length  and 
type  of  firing  line,  and  the  detonator.   For  the  1975  tests  we 
estimated  from  a  comparable  test  setup  that 


FIRING  DELAY  =  +0.31  +  0.05  milliseconds 


(D2) 


The  correction,  VEL.  CORR. ,  to  At     ,  must  be  calculated. 

sound 

It  is  given  by  the  following  equation 


VEL.  CORR.  =  -  --2- 
c0 


r   R/Ro 

u 


L  l 


,R 


Co      /    K, 


_  c°\ 


(D3) 


where 


R0  =  radius  of  the  charge 

c0  =  sound  velocity  in  the  water 

U  =  shock  velocity  in  the  water 

R  =  radial  distance  to  gage 

D  =  detonation  velocity  in  the  charge 

The  first  term  inside  the  brackets  represents  the  contribution  due 
to  the  water  shock  while  the  second,  that  due  to  the  detonation  wave 
Equation  D3  was  derived  from  the  sketch  shown  in  Figure  D-2 . 

The  integrand  in  Equation  D3  was  obtained  by  fitting 

results  from  a  hydrodynamic  code  calculation  for  pentolite  by 

12 
Sternberg  and  Walker     with  the  following  equations 


12/ 


Sternberg,  H.  M. ,  and  W.  A.  Walker,  1971,  "Calculated  Flow  and  Energy 
Distribution  Following  Underwater  Detonation  of  a  Pentolite  Sphere", 
Physics  of  Fluids,  September  1971 
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^-   -1  =  EXP  | 8.84-160.9  +  ^  0%^°  >     for  R/R0  <  25  (D4a) 


—  -1  =  EXP  < -3.264  -  1.098  In  5Z|a.  \  for  r/Rq  >  25*  (D4b) 


To  calculate  VEL.  CORR.  Equation  D3  must  be  integrated  all  the  way 
out  to  the  gage.   (This  correction  does  not  assume  a  constant  value 
beyond  some  arbitrary  point,  i.e.,  for  this  application  the  shock- 
wave  cannot  be  considered  to  be  like  a  sound  wave,  anywhere.) 

Table  D-l  summarizes  the  results  of  the  sound  ranging 

computations.   The  last  two  columns  give  the  comparisons  for  the 

time  of  arrival  for  the  direct  Shockwave.   The  comparison  is  easily 

understood  by  subtracting  At     ,  from  both  sides  of  Equation  Dl . 
2  ^         sound 

The  maximum  discrepancy  is  0.4  milliseconds  (Shot  786)  or  a  spatial 
discrepancy  of  about  2  feet. 

Tables  D-2  thru  D-7  list  cage  and  gage  locations  for  the 
19  75  test  series  which  where  calculated  by  the  method  of  this  appendix, 
The  values  listed  in  Tables  D-2  thru  D-7  were  used  for  the  present 
analysis  of  the  1973  and  1975  test  series  data.   The  distance  along 
the  support  wire  (Column  1)  is  the  nominal  or  desired  depth  of  the 
fish  cage  or  pressure  gage.   The  wire  angle  (Column  2)  is  the 
calculated  deviation  from  the  vertical  of  the  wire  segment  attached 
to  the  preceding  depth  coordinate.   The  cage  and  gage  coordinates 
(Columns  3  thru  6)  are  calculated  relative  to  the  point  on  the  water 
surface  directly  above  the  charge.   Sound  speeds  were  calculated 

from  salinity  and  water  temperature  at  charge  depth  using  Del  Grosso's 

13 
equation.    (All  the  computations  described  in  this  appendix  were 

done  using  an  HP-65  Programmable  Pocket  Calculator.) 


*This  straight  line  log-log  extrapolation  beyond  the  range  of  Sternberg  and  Walker's 
calculation  should  be  adequate  out  to  at  least  R/R0  =  1000. 

13. Del  Grosso,  V.  A.,  "New  Equation  for  the  Speed  of  Sound  in  Natural  Waters  (with 
Comparisons  to  Other  Equations)",  J.  Acoust .  Soc.  Am.,  October  197*+ 
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TABLE  D-2 
CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  782 

Charge:  70.4   lbs  pentolite 
Charge  wire  length:  30  feet 
Charge  wire  angle:   0  degrees 
Charge  depth:   30  feet 

Charge-to-gage  line  distance:  300  feet 
Sound  speed:   48  50   ft/sec 
Water  depth*:   156  feet 


WIRE 
DISTANCE 

WIRE 
ANGLE 

CAGE  AND 

GAGE  COORDINATE 

S 

HOR. 

RANGE 

DEPTH 

HOR. 

RANGE 

DEPTH 

(FEET) 

(DEG.) 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 



300 

0.0 

91.4 

0.0 

5 

0 

300 

5.0 

91.4 

1.5 

10 

0 

300 

10.0 

91.4 

3.0 

15 

0 

300 

15.0 

91.4 

4.6 

40 

0 

300 

40.0 

91.4 

12.2 

45 

0 

300 

45.0 

91.4 

13.7 

50 

0 

300 

50.0 

91.4 

15.2 

55 

0 

300 

55.0 

91.4 

16.8 

57.5 

0 

300 

57.5 

91.4 

17.5 

77.5 
87.5 

0 
0 

300 
300 

77.5 
87.5 

91.4 
91.4 

23.6 
26.7 

97.5 

0 

300 

97.5 

91.4 

29.7 

♦Calculated  from  times  of  arrival  of  bottom  reflection 
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TABLE  D  -3 
CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  783 

Charge:    70.2  lbs  pentolite 
Charge  wire  length:  30   feet 
Charge  wire  angle:    0  degrees 
Charge  depth:  30   feet 

Charge-to-gage  line  distance:  200   feet 
Sound  speed:  4857    ft/sec 
Water  depth*:   152  feet 


WIRE 
DISTANCE 

WIRE 
ANGLE 

CAGE  ANE 

)  GAGE  COORDINATES 

HOR. 

HOR. 

RANGE 

DEPTH 

RANGE 

DEPTH 

(FEET) 

(DEG.) 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 



200.0 

0.0 

61.0 

0.0 

5 

-10.5 

199.1 

4.9 

60.7 

1.5 

10 

-10.5 

198.2 

9.8 

60.4 

3.0 

15 

-10.5 

197.3 

14.7 

60.1 

4.5 

40 

-10.5 

192.7 

39.3 

58.7 

12.0 

45 

-10.5 

191.8 

44.2 

58.5 

13.5 

50 

-10.5 

190.9 

49.2 

58.2 

15.0 

55 

-10.5 

190.0 

54.1 

57.9 

16.5 

60 

-10.5 

189.1 

59.0 

57.6 

18.0 

80 

-10.5 

185.4 

78.7 

56.5 

24.0 

90 

-  8 

184.0 

88.6 

56.1 

27.0 

100 

-  5 

183.2 

98.5 

55.8 

30.0 

*Calculated  from  times  of  arrival  of  bottom  reflection. 
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TABLE    D-4 

CORRECTED  CAGE  AND  GAGE  LOCATIONS 

SHOT   78  4 

Charge:   71.6   lbs  pentolite 
Charge  wire  length:  30   feet 
Charge  wire  angle:    0  degrees 
Charge  depth:  30   feet 

Charge-to-gage  line  distance:  300   feet 
Sound  speed:  4866    ft/sec 
Water  depth*:  135   feet 


WIRE 
DISTANCE 

WIRE 
ANGLE 

CAGE  ANC 

GAGE  COORDINATES 

HOR. 

HOR. 

RANGE 

DEPTH 

RANGE 

DEPTH 

(FEET) 

0 

(DEG.) 

(FEET) 
300.0 

(FEET) 
0.0 

(METERS) 
91.4 

(METERS) 
0.0 

5 

11 

301.0 

4.9 

91.7 

1.5 

10 

11 

301.9 

9.8 

92.0 

3.0 

20 

10 

303.6 

19.7 

92.5 

6.0 

30 

10 

305.4 

29.5 

93.1 

9.0 

40 

9 

306.9 

39.4 

93.5 

12.0 

45 

8 

307.6 

44.3 

93.8 

13.5 

50 

6 

308.2 

49.3 

93.9 

15.0 

55 

4 

308.5 

54.3 

94.0 

16.6 

60 

3 

308.8 

59.3 

94.1 

18.1 

80 
90 

2 

1 

309.5 
309.6 

79.3 
89.3 

94.3 
94.4 

24.2 
27.2 

100 

1 

309.8 

99.3 

94.4 

30.3 

♦Calculated  from  times  of  arrival  of  bottom  reflection. 
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TABLE  D-5 
CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  785 

Charge:  1.25    lbs  pentolite 
Charge  wire  length:  30   feet 
Charge  wire  angle:  -5   degrees 
Charge  depth: 29.9  feet 

Charge-to-gage  line  distance:  40    feet 
Sound  speed:  4864    ft/sec 
Water  depth*:  160   feet 


WIRE 
DISTANCE 

WIRE 
ANGLE 

CAGE  ANE 

)  GAGE  COORDINATES 

HOR. 

HOR. 

RANGE 

DEPTH 

RANGE 

DEPTH 

(FEET) 

(DEG.) 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 



42.6 

0.0 

13.0 

0.0 

5 

-11 

41.7 

4.9 

12.7 

1.5 

10 

-11 

40.7 

9.8 

12.4 

3.0 

20 

-11 

39.3 

19.7 

12.0 

6.0 

30 

-  7 

38.1 

29.6 

11.6 

9.0 

40 

-  7 

37.1 

39.6 

11.3 

12.1 

50 

-  5 

36.0 

49.5 

11.0 

15.1 

60 

-  4 

35.0 

59.5 

10.7 

18.1 

70 

-  2 

33.9 

69.4 

10.3 

21.2 

80 

+  2 

34.1 

79.4 

10.4 

24.2 

90 

+  5 

34.4 

89.4 

10.5 

27.2 

100 

+  5 

35.1 

99.4 

10.7 

30.3 

*Calculated  from  times  of  arrival  of  bottom  reflection, 
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TABLE    D-6 

CORRECTED  CAGE  AND  GAGE  LOCATIONS 

SHOT   786 

Charge:   1.26  lbs  pentolite 
Charge  wire  length:   30  feet 
Charge  wire  angle:  -10  degrees 
Charge  depth:  29. 5  feet 

Charge-to-gage  line  distance:   40  feet 
Sound  speed:    4886  ft/sec 
Water  depth*:   115  feet 


WIRE 
DISTANCE 

WIRE 
ANGLE 

CAGE  ANC 

GAGE  COORDINATES 

HOR. 

HOR. 

RANGE 

DEPTH 

RANGE 

DEPTH 

(FEET) 
0 

(DEG.) 

(FEET) 
45.2 

(FEET) 
0.0 

(METERS) 
13.8 

(METERS) 
0.0 

5 

-30 

42.7 

4.3 

13.0 

1.3 

10 

-30 

40.2 

8.7 

12.3 

2.7 

20 

-24 

36.1 

17.8 

11.0 

5.4 

30 

-19 

32.9 

27.3 

10.0 

8.3 

40 

-15 

30.3 

36.9 

9.2 

11.2 

50 

-12 

28.2 

46.7 

8.6 

14.2 

55 

-10 

27.4 

51.6 

8.4 

15.7 

60 

-10 

26.5 

56.5 

8.1 

17.2 

65 

-  8 

25.8 

61.5 

7.9 

18.7 

70 

-  6 

25.3 

66.5 

7.7 

20.3 

75 

-  4 

24.9 

71.5 

7.6 

21.8 

80 

-  2 

24.7 

76.4 

7.5 

23.3 

♦Calculated  from  times  of  arrival  of  bottom  reflection. 
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TABLE  D  -7 
CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  787 

Charge:   72.2   lbs  pentolite 
Charge  wire  length:   10  feet 
Charge  wire  angle:   0   degrees 
Charge  depth:   10  feet 

Charge-to-gage  line  distance:   30  0  feet 
Sound  speed:  4922    ft/sec 
Water  depth*:  14  8   feet 


WIRE 
DISTANCE 

WIRE 
ANGLE 

CAGE  ANE 

i  GAGE  COORDINATES 

HOR. 

HOR. 

RANGE 

DEPTH 

RANGE 

DEPTH 

(FEET) 

(DEG.) 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 



300 

0 

91.4 

0.0 

5 

0 

300 

5 

91.4 

1.5 

10 

0 

300 

10 

91.4 

3.0 

40 

0 

300 

40 

91.4 

12.2 

50 

0 

300 

50 

91.4 

15.2 

55 

0 

300 

55 

91.4 

16.8 

60 

0 

300 

60 

91.4 

18.3 

70 

0 

300 

70 

91.4 

21.3 

80 

0 

300 

80 

91.4 

24.4 

90 

0 

300 

90 

91.4 

27.4 

100 

0 

300 

100 

91.4 

30.5 

*Calculated  from  times  of  arrival  of  bottom  reflection 
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APPENDIX  E 
NOTE  ON  FISH  CLOSE  TO  THE  WATER  SURFACE 


According  to  the  bladder  oscillation  model  for  explosion 
injury,  as  the  fishes'  depth  tends  to  zero  near  the  water  surface 
the  kill  probability  should  also  tend  to  zero  reqardless  of  how  high 
the  incident  pressure.   Eventually,  however,  for  fish  at  the  water 
surface  and  sufficiently  close  to  the  charge,  the  bladder  oscillation 
model  must  break  down  and  a  different  mechanism  for  injury  must  take 
over.   Apparently,  this  happens  so  close  to  the  charge,  and  to  the 
water  surface,  to  be  of  little  practical  importance,  since  results 
reported  by  the  Lovelace  Foundation    show  no  evidence  of  such  a 
transition  for  149  gm  Carp  (length  =  21  cm)  at  a  depth  of  5  cm 
subjected  to  5.58  megapascals  (810  psi)  incident  pressure. 

A  possible  close-in  non-bladder  injurv  mechanism  is  tissue 
damage  due  directly  to  compression  by  the  incident  Shockwave.   If 
this  is  the  case,  tissue  damage  probably  also  occurs  due  to  the 
pressure  wave  (generally  an  order  of  magnitude  areater)  emitted  by 
the  oscillating  swim  bladder.   Assuming  that  this  is  the  meaning  of 
the  ratio,  AMIN/A.,  in  the  damage  parameter  Z,  we  rewrite  Equation 
3.1.3  in  terms  of  the  adiabatic  compression  ratio,  PMIN/p0 ,  as 


__1 
/pmtn\ 

X 


=  _100  m  (^_N)  3Y  (ei) 


where  PMIN  is  the  maximum  value  of  the  oscillating  bladder  pressure 
(corresponding  to  radius  AMIN  at  the  first  compression) ,  p^  is  the 
initial  ambient  pressure,  and  y, the  adiabatic  exponent  for  air. 
Using  El  we  then  rewrite  3.1.5  as 

E-l 
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Z  =  i°-°  In  PMIN/^  +  100  In  AMAX/A±  (E2) 


Note  that  Equation  E2  is  equivalent  to  3.1.5.   It  is  just  our 
point  of  view  that  has  changed. 

Now  let  us  suppose,  for  fish  very  close* to  the  water 
surface,  that  the  Shockwave  peak  overpressure  PMAX  is  the  only- 
tissue  damaging  pressure.   The  damage  parameter  by  this  mechanism 
equivalent  to  Equation  E2  is  then  given  by 


100  , 
PMAX  '  3y 


^+1  (E3) 

.  Pi      J 


where  PMIN  is  replaced  by  PMAX  +  p.  and  the  second  term  drops  out. 

Note  that  Z _,..,..,  must  always  be  smaller  than  Z  calculated  by 
PMAX  J 

Equation  3.1.5  except  for  those  cases — the  ones  of  interest  here-- 

where  the  pressure  wave  is  of  such  short  duration  that  the  transient 

response  of  the  swimbladder  is  suppressed*.   Thus,  when  in  doubt 

which  damage  parameter  to  use — zPMAv  (Equation  E3)  or 

Z  (Equation  3 . 1 . 5) --calculate  both  parameters,  and  use  the  greater 

value. 

Finally,  we  estimate  the  overpressure  level,  PMAX, 
associated  with  the  assumption  that  all  the  damage  corresponding 
to  the  damage  parameter  Z     is  due  directly  to  tissue  compression 
by  the  incident  Shockwave.   Solving  E3  for  PMAX  we  get 

"* 

Since  Z       is  equivalent  to  the  damage  parameter  X 

(Equation  3.1.3)  calculated  from  the  non-oscillatory  steady-state 
response  of  an  infinitely  small  bubble  to  an  overpressure  PMAX 
of  finite  rise  time. 

E-2 
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PMAX  -  Pi    EXP[  (3yZpMAX/100)  -  1]  (E4) 


Setting  p.  =  1  atmosphere  (1.014  x  10   pascals)  and 

Z™,,^,  =  125  (50%  kill  value—Equation  3.2.1)  we  get 
PMAX 


PMAX  =19.2  MPa  (=2790  psi)  (E5) 

(approximately  24  charge  diameters) 

for  the  overpressure  level  PMAX  corresponding  to  50%  kill 
by  this  mechanism  to  fish  which  are  very  close  to  the  water 
surface . 


E-3/E-4 
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